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编译 者 序 


“7 沁 本 书 主要 介绍 电子 计算 机 解 各 类 方程 的 方法 。 全 忆 


坟 逐 次 逼近 为 中 心 ， 系 统 地 论述 了 高 次 方程 、 无 理 方 
程 、 线 性 方程 组 、 非 线性 方程 组 等 应 用 问题 的 计算 机 解 
法 ， 对 每 种 方法 都 列 出 了 具体 求解 公式 ， 并 作 了 直观 的 
几何 解释 ， 最 后 还 简单 讨论 了 方法 的 妆 效 性 。 

本 书 的 特点 是 主题 突出 、 邮 材 精 当 ， 与 计算 机 解 题 
联系 紧密 ， 着 重 方 法 的 实用 性 ， 完 全 还 开 了 用 高 等 数学 
知识 进行 此 琐 复 杂 的 理论 推导 ， 而 从 说 茜 数学 出 发 ， 形 
象 直观 地 对 各 种 方法 展开 讨论 。 会 书 以 通俗 的 文字 ， 从 。 
”生活 中 实际 例子 出 发 ， 由 浅 入 深 地 未 步 把 用 计算 机 解 题 
的 各 种 方法 呈现 给 读者 。 因 此 ， 具 有 有 一般 中 学 数学 基础 . 
如 识 的 人 部 能 阅读 这 本 书 ， 并 能 掌握 书 中 论述 的 方法 。 

本 书 系 根据 莫斯科 罗 蒙 诺 索 夫 国立 大 学 的 N*Ya* 维 
伦 金 先生 的 同名 著作 编译 而 成 。 原 著 用 依 文 出 版 过 两 次 ， 
后 被 译 成 英文 ， 是 一 本 很 有 特色 的 计算 机 应 用 数学 的 曾 
及 性 读物 。 可 作为 学 习 计 算 机 解 题 方 法 的 启蒙 教材 ， 也 
可 供 从 事 计 算 机 工作 的 工程 技术 人 员 参 考 。 还 可 作 数 学 
爱好 者 的 参考 读物 。 

在 本 书 编译 过 程 中 ， 斑 仲 文 先生 审 校 了 全 部 稿件 并 
对 文章 作 了 不 少 修改 ; 张 新 副 教授 也 审 闵 了 全 部 文稿 。 
对 他 们 给 予 的 极 大 的 帮助 和 支持 ， 作 者 庶 在 此 致谢 。 邮 
于 作者 水 平 有 限 夫 译 误 之 处 难免 ， 攻 请 读者 指正 。 


潘 林 森 王 飞 
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1. 绪 言 


在 中 学 学 习 数 学 的 时 候 ， 我 们 大 量 的 时 间 是 花 在 求解 方 
程 或 方程 组 上 ， 最 初学 习 的 是 一 次 方程 和 一 次 方程 组 丑 
来 ， 又 学 习 了 二 次 方程 、 双 二 次 方程 和 无 理 方程 ， 最 后 ， 通 ， 
熟悉 和 掌握 了 指数 方程 、 对 数 方程 和 三 角 方 程 。 & 

我 们 如 此 重视 方程 决 不 是 偶然 的 ， 因 为 方程 在 数学 实际 
应 用 中 是 非常 重要 的 。 在 你 所 选择 的 任何 一 个 领域 中 应 用 数 
学 ， 都 要 求解 方 垂 或 方程 组 ， 才 能 得 到 最 终 的 答案 。 

在 中 学 ， 方 程 常 常用 于 求解 物理 问题 ， 例 如 ， 考 虑 下 面 
的 问题 ; 

把 一 石头 扔 进 井 轧 ， 如 果 石 头 扔 下 7? 秒 后 便 听 到 石头 磁 
击 水 面 的 声音 ， 求 并 水 水 面 至 井 日 的 距离 ? 

我 们 假设 距离 为 x， 为 此 得 到 一 个 关于 x 的 方程 


/2+2=7 
R y 
其 中 ，V 是 声音 在 空中 传播 的 速度 ( 2578 是 石头 洲 到 间 


水 水 面 的 时 间 ，x/v 是 声音 从 井 水 水 面 传 到 地 面 的 时 间 )。 
这 是 一 个 无 理 方程 。 令 wx 二 y， 原 式 化 为 一 个 二 次 方程 


yg /2 +2- 
+ 2 7 =0 


这 个 方程 ， 用 大 家 部 知道 的 二 次 方程 求解 公式 便 可 求解 。 
Pd 和 
1 


方程 同样 可 用 于 解答 风 何 问题 ， 例 如 把 一 长 度 为 ! 的 线 
段 4B 划 分 成 4C 和 CB， 使 得 4B : 4C=4C : CB， 这 个 间 
题 便 导 出 一 个 二 次 方程 

x*+lix—1:= 0 
其 中 ，x 表 示 AC 的 长 度 。 

把 一 角度 为 a 的 角 分 为 3 等 份 的 问题 又 导 出 一 个 比较 复 

杂 的 方程 , 


42Xs 一 3x 一 cosa 一 0 
其 中 ，x 一 cos3。 这 样 的 方程 时 三 次 方程 。 在 中 学 是 没有 


学 习 过 的 ， 科大 都 证 明了 关于 三 
， 次 方程 有 求解 公式 ， 见 下 面 的 公式 ( 3 ) 。 

然而 ， 在 物理 中 ， 我 们 常常 遇 到 的 问题 ， 导 出 的 是 更 复 
泰 的 方程 ， 这 些 方程 ， 往 往 用 我 们 过 去 在 中 学 或 大 学 所 学 的 
知识 是 难以 求解 的 ， 例 如 ， 把 一 根 铁 条 (工程师 称 它 为 雁 
! 条 ) ， 铅 直 地 固定 紧 它 的 一 端 ， 这 时 ， 如 果 我 们 敲 击 铁 条 一 
下 ， 铁 条 就 要 产生 横 波 振动 ， 要 求 这 种 振动 的 频率 ， 就 要 求 
解 方程 


3 ~ 2 
ste cosz (1) 


其 中 ，e 二 2.71828…… 

中 学 里 ， 根 本 没有 求解 这 种 方程 的 法 则 ， 不 要 以 为 这 是 
由 于 中 学 数学 课程 太 简短 ， 实 际 上 ， 在 中 学 通常 所 学 的 知 
识 范 围 内 ,方程 ( 1 ) 是 根本 没有 公式 求解 的 。 现 在 ， 让 我 
们 对 这 个 问题 作 比 较 确 切 的 陈述 : 

一 个 方程 ， 如 果 它 的 根 能 够 借助 于 算术 运算 、 开 方 运算 
自 及 指数 、 对 数 、 三 角 函 数 、 反 三 角 函 数 ， 用 方程 的 各 种 参 
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数 来 表示 。 我 们 就 说 方程 有 求解 公式 。 在 这 种 意义 上 ， 二 次 
方程 22 十 xP 十 4= 0 有 以 下 求解 公式 


x 一 -一 去 二 v374 一 9 (2) 
桐 样 ， 三 次 方程 " 
X 1 十 4 一 0 
也 有 求解 公式 
下 _ 了， /二 -下 一 全- gp 
= Lt ++ '/-! EE + 好 (3) 


可 是 在 实际 使 用 公式 ( 3 ) 时 ， 却 包含 了 若干 困难 ， 同 时 通 
要 求 使 用 复数 。 
对 于 四 次 方程 ， 也 有 求解 公式 ， 但 非常 复杂 ， 因 此 ， 赛 
们 不 在 此 讨论 7。 
五 次 或 高 于 五 次 方程 的 情形 就 更 糖 了 ， 挪 威 数学 家 阿 员 
尔 早 在 1826 年 就 证 明了 : n 次 代数 方程 : 
doXn 十 01X2 1 十 wn 十 an 一 0 


当 m 之 5 时 ， 没 有 求解 公式 〈 它 的 根 不 能 借助 于 算术 运算 和 
了 研 方 运算 来 表示 ) ， 只 是 在 特殊 情况 下 ， 高 于 四 次 的 代数 
方程 才 可 能 有 求解 公式 “。 

如 果 数 学 家 们 研究 方程 时 仅仅 限于 考虑 准确 解 ， 即 只 考 


。 任何 三 次 方程 8ox3 十 alx2 十 az2x 十 a3= 0 都 可 借助 于 代 换 x 十 5 
Al/350=Y 化 成 上 而 的 形式 。 让 


。 关于 代数 方程 可 参见 A.G. Kurosh, Algebraic Equstions of 
arbitrary Degreces，MirPublishers，Moscow，1977 


虑 能 由 求解 公式 表示 的 解 ， 那 么 ， 工 程 师 和 数学 家 之 间 将 会 
发 生 下 面 的 谈话 : 

工程 师 ， 我 设计 一 村 建筑 时 ， 遇 到 这 个 方程 ( 拿 出 方 
程 ) 需要 尽快 地 知道 它 的 解 一 一 因为 一 个 月 之 内 我 必须 完成 
任务 。 

数学 家 ， 我 很 高 兴 帮 助 你 ， 不 过 ， 这 种 形式 的 方程 没有 
求解 方法 。 

工程 师 ， 你 不 能 推导 一 个 求解 公式 吗 ? 

数学 家 ， 那 可 没有 用 处 ， 很 久 以 前 就 证 明了 这 种 方程 没 
有 求解 公式 。 

昕 了 这 样 的 谈话 ， 人 们 可 以 想像 ， 工 程 师 的 数学 见解 将 
使 方程 求解 的 可 能 人 性 交 得 更 粮 。 幸 运 的 是 ， 不 可 能 发 生 这 样 
的 谈话 。 实 际 上 ， 工 程 师 通常 不 需要 这 样 那 样 方程 的 求解 公 
式 。 他 们 所 需要 的 是 满足 一 定 精 度 的 正确 答案 一 一 而 这 种 答 
案 或 者 用 公式 求 得 或 者 用 别 的 方法 求 得 对 他 们 者 不 太 感 兴 
趣 。 

比如 ， 我 们 设想 由 某 个 公式 得 出 的 答案 是 xz 一 3 十 Vi5。 
很 明显 ， 这 个 答案 不 能 直接 用 在 实际 中 ( 人 们 不 可 能 在 机 器 
上 制造 出 精确 长 度 为 ( 3 十 13 ) cm 的 零件 ) ， 实 际 应 用 
中 ，v/ 王 应 该 表示 成 十 进 制 小 数 ， 小 数 的 位 数 要 由 给 出 的 实 
际 问题 的 要 求 精度 而 决定 。 

因此 ， 只 要 数学 家 告诉 工程 师 怎样 计算 方程 的 根 以 及 必 
需 的 精度 ， 工 程 师 就 非常 满意 了 。 激 学 中 已 呈现 出 了 大 量 的 
求 方程 近似 解 的 方法 ， 本 书 中 将 要 氢 述 其 中 一 些 方法 。 


2. 逐 次 遥 近 


大 多 数 求 方程 近似 解 的 方法 都 是 报 据 逐次 逼近 的 概念 ， 
这 种 概念 不 仅 可 用 于 求解 方程 ， 而 且 也 可 以 用 来 解决 许多 实 
际 问 题 。 

逐次 逼近 法 或 尝试 法 及 误差 法 常常 被 炮 手 们 采用 。 他 们 
为 了 击 中 车子， 需要 确定 大 炮 射 击 的 方位 和 目标 ， 炮 弹 发 出 
居 ， 观 察 炮弹 的 爆炸 点 ， 如 果 没 击 中 靶子 ， 就 适当 修正 大 炮 
射击 的 方位 和 目标 ， 然 后 再 继续 射击 。 经 过 几 次 这 样 调 整 逼 
近 以 后 ， 他 们 就 能 确定 出 准确 的 方位 和 目标 ， 以 便 击 中 葛 
子 。 


图 1 


有 时 候 ， 决 定 瞻 准点 也 需要 用 逐次 逼近 法 ， 比 如 ， 一 加 
离 射 炮 在 O 点 射击 一 架 飞 机 ( 图 1 ) 。 假 定 某 一 时 刻 ， 飞 机 
在 4 点 ， 此 时 如 果 高 射 雹 向 4o 点 瞄准 射击 ， 是 击 不 中 .飞机 
的 。 因 为 炮弹 射击 的 同时 ， 飞 机 已 飞 到 了 41 点 。 如 果 大 位 
知道 了 飞机 的 速度 和 炮弹 的 速度 ， 就 可 比较 容易 找到 射击 此 
准点 41!,， 可 是 ， 如 果 仍 然 瞻 准 4 点 射击 ， 还 是 击 不 中 飞机 
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的 ， 因 为 炮 简 的 倾斜 度 改变 了 炮弹 运行 的 轨道 ， 因 此 ， 炮 弹 
在 04。 距 离间 运行 的 时 间 与 在 04, 距 离间 运行 的 时 间 不 一 
样 ， 结 果 ， 炮 弹 伊 击 不 中 飞机 。 但 是 ， 了 瞄准 A: 点 射击 就 比 
足 准 4 点 射击 失误 小 。 为 了 使 失 变 的 可 能 性 更 小 ， 炮 手 们 
还 应 该 求 得 炮弹 在 OA 距离 间 运 行 的 时 间 以 及 在 这 个 时 间 内 
飞机 所 到 达 的 地 点 。4: 点 将 是 跟踪 目 标点 的 再 次 逼 近 ， 之 
后 再 根据 炮弹 由 O 点 到 4; 点 所 需 时 间 来 计算 4 点 ， 即 是 飞 
机 在 那个 时 间 内 要 到 达 之 处 ， 经 过 车 干 次 允 近 以 后 ， 我 们 将 
找到 具有 一 定 准 确 度 的 瞄准 点 。 

逐次 逼近 的 方法 还 可 用 来 解决 许多 其 它 问 题 。 

假定 ， 人 们 要 从 N 个 采 沙 场 4,、4,…… ，4, 运 沙 到 于 
个 建筑 工地 Bi， Bs, ee ， Bn, 假设 第 /个 沙场 4: 每 天 生 
产 aj 吨 沙 ， 第 K 个 工地 B4 每 天 需要 尼 吨 沙 ， 再 假设 从 4y 沙 
场 运输 一 吨 沙 到 B4 工 地 需要 的 费用 为 Cyr( 由 4, 到 Bs 之 间 的 
距离 和 运输 单价 等 决定 ) 。 

为 了 制定 运输 计划 ， 我 们 作 表 1 ， 在 表 1 中 ，xy: 表 示 
从 4 ;沙场 运 到 Bx 工 地 沙 的 总 数 。 


表 1 

| Bi Bs | Be | mm 

A X11 X12 1k Xim 
A, | P21 ; | 


黄 中 ，*jx 满 足 关系 : 


Xi 十 Xi 十 十 zjim 和 ai 
《每 天 从 4 ;沙场 运 出 的 沙 石 不 能 超过 4a; 号 ,二 ，. 吉 : 
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Kk Xr 二 十 xs 一 中 


(《 Bs 工地 每 天 要 得 到 纪 吨 沙 ) 
如 果 按 表 1 的 运输 计划 ， 其 运输 费用 为 : 
C 一 Cuixii 十 Cisxis 十 ，…… 十 Caxkte 十 
十 Cliyx:i 十 Caxs: 十 和 十 Csxx。 十 
十 Caixal 十 Cnsxms 十 …… 十 Cmaxms (4》 


任何 运输 计划 都 应 该 是 使 费用 C 尽 可 能 最 少 ， 首 先 ， 
我 们 试探 性 地 拟定 一 种 计划 ， 比 如 ， 可 采用 下 面 的 运输 方 
案 : 

开始 ， 把 4, 沙 场 的 沙 运 到 离 它 最 近 的 工地 ,如 果 这 个 沙 
场 每 天 生产 的 沙 除了 供应 那个 工地 外 还 有 剩余 ， 就 把 剩 下 的 
沙 运 到 其 它 工地 中 离 它 最 近 的 一 个 工地 上 ,如 还 有 剩 ,再 把 简 
下 的 运 到 离 它 较 近 另 一 个 工地 。4, 沙 场 生产 的 沙 运 完 以 后 ， 
再 同样 安排 运 4, 沙场 的 沙 ， 如 此 等 等 。 直 到 最 后 ， 每 个 工 
地 都 有 限定 的 沙场 运 沙 。 

可 是 ， 这 种 运输 计划 并 不 是 最 好 的 方法 ， 因 为 这 样 安排 
运 沙 ， 可 能 剩 下 最 后 几 个 建筑 工地 要 到 很 远 的 沙场 去 运 沙 ，。 
所 以 ， 这 种 运输 计划 必须 加 以 改进 。 有 些 人 采用 的 运输 方案 
是 :在 沙场 和 工地 之 语 ， 先 考虑 给 需 沙 量 大 的 工地 运输 ， 再 
考虑 给 需 沙 量 小 的 工地 运输 。 

改变 计划 使 得 费用 减少 ， 这 种 方法 是 数学 的 一 个 分 支 ， 
则 线性 规划 *。 

借助 于 这 些 方法 ， 经 过 几 次 逼近 后 ， 我 们 就 能 找 出 一 种 


?关于 线性 规划 可 人 参见 A.S.Solodovnikov 
Introductoni to linear algebra and linear programming” ,. 


Prosveshcheniye, Moscow, 1966 


分 


使 和 式 ( 4 ) 的 值 最 小 或 者 接近 最 小 的 计划 。 

通常 ， 当 我 们 制定 一 种 计划 、 一 张 时间 硝 等 等 的 时 候 ， 
最 初 总 是 用 一 些 比较 粗略 的 近似 值 ， 跑 后 再 逐次 改进 ， 直 到 
得 出 最 终 的 结果 为 止 。 

在 工厂 车 间 里 ， 生 产 加 工 机 器 零件 ， 也 可 以 署 成 是 一 种 
逐次 腺 近 的 过 程 。 开 始 是 一 种 粗 覆 的 逼近 一 一 拿 一 块 铸件 或 
晴 件 ， 再 把 坯 件 拿 到 车 床上 先 制 成 近似 零件 形状 ， 然 后 在 更 
精密 的 车 床上 经 过 几 次 加 工 ， 也 就 是 几 次 逼近 以 后 ， 便 生产 
出 我 们 需要 的 机 器 零件 。 


3.. 阿 基 里 斯 和 乌龟 


最 先 提 到 逐次 逼近 法 的 是 公元 500 年 前 伊 利 亚 的 地 诺 。 
这 位 哲学 家 试图 证 明 自 然 界 没有 运动 ， 他 用 下 面 的 论据 证 明 
不 存在 运动 ， 即使 让 跑 得 最 癸 的 希腊 神 阿 基 里 斯 去 追赶 工 只 
乌龟， 他 也 追 不 上 。 先 假设 阿 基 里 斯 和 乌龟 之 间 相 距 1000 
步 远 距离 ， 和 再 假设 阿 基 里 斯 1 秒 钟 跑 10 步 ,乌龟 1 秒 钟 谍 动 
1 步 ， 那 么 ， 在 100 秒 钟 时 ， 阿 基 里 斯 跑 完 他 与 乌龟 相距 的 
1000 步 距离 ， 但 与 此 同时 ， 鸟 急 也 胞 了 100 步 远 。 再 用 10 种 
钟 时 间 ， 阿 基 里 斯 又 路 完 这 100 步 远 的 距离 ,但 乌龟 同时 也 又 
爬 了 10 步 远 ， 等 阿 基 轩 斯 再 哆 完 这 10 步 时 ,乌龟 又 私 了 工 步 。 
因此 ， 乌 兔 永远 在 阿 基 轩 斯 前 面 ， 阿 基 里 斯 永远 不 能 赶 上 岛 
龟 。 所 以 ， 自 然 界 没有 运动 。 

当然 ， 季 诺 的 论证 仅仅 是 一 种 措 词 巧妙 的 雇 论 ， 是 没有 
机 由 的 ， 运 动 本 身 是 物体 内 在 的 性 质 。 

要 计算 阿 基 硬 斯 追赶 上 乌龟 的 时 间 ， 对 任何 小 学 生 都 不 
是 困难 的 事 ， 建 立 一 个 方程 式 就 可 计算 ， 

10x 一 zx 一 1000 (5) 

其 中 ，x 是 所 求 时 间 ， 由 方程 我 们 得 到 


不 过 ， 季 诺 的 论证 可 以 看 成 是 一 种 近似 求解 方程 ( 5 》 
的 特殊 方法 。 
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让 


实际 上 ， 把 x 项 移 到 方程 的 右边 ， 两 边 再 除 以 10， 我 们 
将 得 到 方程 
100 二 < 
*=100+ 人 (6) 


和 如果 有 略 去 方程 右边 的 x*/10 项 ( 它 相 对 于 x 是 小 量 ) ， 我 们 得 
到 x 的 近似 解 2: 二 100。 现 在 ， 把 近似 解 *, = 二 100 代 入 方程 右 
边 的 x*， 我 们 将 得 到 一 个 更 精确 的 x* 值 ， 即 x, = 100 十 10 一 
110， 再 把 这 个 新 值 代 入 方程 右边 *， 又 得 到 下 一 个 近 似 解 
za 一 100 十 110V10 一 111， 按 这 种 方法 ， 我 们 得 到 一 组 近似 
解 
%1 一 100，X2 一 110，%s 一 111，X4 一 111.1， 

也 就 是 说 ， 季 诺 的 论证 得 出 了 上 述 这 组 数 。 这 组 数 能 够 通过 
关系 式 


xati=100+xs/10 (7) 

对 "逐次 计算 得 出 。 当 x* 增 大 时 ， 这 组 数 逼 近 方程 ( 5 ) 的 准 
=111 工 

确 解 *=1l1g 


因为 x/10 项 较 之 x* 是 小 量 ， 所 以 上 面 描 述 的 求解 方法 是 
成 功 的 。 不 然 的 话 ， 我 们 将 得 到 一 组 越 来 越 离 准确 解 远 的 
数 。 例 如 ， 假 设 阿 基 里 斯 不 是 追赶 一 只 息 入 很 慢 的 乌 旬 ， 而 
是 追赶 一 只 跑 得 很 快 的 羚羊 ， 我 们 假设 凑 羊 每 秒 钟 跑 20 步 
远 ， 为 求 阿 基 里 斯 追 上 冷 羊 所 需 的 时 间 ， 我 们 必须 解 方程 
10x 一 20x 一 1000 (8) 
这 个 方程 的 准确 解 是 x*= 一 100。 这 就 是 说 ， 阿 基 里 斯 和 券 单 
并 驾 齐 上 驱 的 时 间 是 在 100 秘 钟 以 前 ， 现 在 券 羊 已 经 跑 在 阿 雅 
里 斯 的 前 面 了 ， 它 们 相隔 的 距离 随 着 时 间 增 加 在 增 大 。 
现在 ， 我 们 试 着 把 解 方程 ( 5 ) 的 方法 用 来 解 方程 (8 )， 
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与 解 方程 ( 5 ) 一 样 ， 先 把 20x 移 到 方程 右边 ， 然 后 ， 两 边 
周 除 以 10， 我 们 得 到 方程 

X 一 100 十 2x (9 ) 
在 方程 右边 令 *, 二 0， 得 到 x, 二 100, 将 这 个 值 代 入 方程 (9 ) 
右边 得 到 第 二 个 近似 解 x*, 二 300， 继 续 这 个 过 程 ， 便 得 到 了 
一 组 数 

Xo=0, xX!1=100, 7#, 一 300，Xs 一 700， 
我 们 看 到 ， 这 组 数 是 不 逼近 于 方程 ( 8 ) 的 准确 解 z 一 一 100 
的 。 
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4 ， 计 算 机 做 除法 


读者 可 能 感到 迷惑 的 是 ， 上 既然 求解 方程 ( 5 ) 的 准确 解 
沉 很 简单 的 ,为 什么 我 们 还 要 用 逐次 逼近 法 去 求解 呢 ? 当 然 
方程 ( 5 ) 本 身 对 我 们 自然 没有 多 大 兴趣 。 我 们 感 兴趣 的 是 
这 种 逐次 逼近 方法 ， 目 的 是 要 把 这 种 方法 应 用 到 更 复杂 的 方 
程 上 去 。 

值得 指出 的 是 ， 随 着 高 速 电 子 计算 机 的 出 现 ， 用 逐次 到 
近 法 求解 与 ( 5 ) 相 类 似 的 方程 ， 目 前 变 得 很 必要 了 。 有 此 
计算 机 只 能 进行 三 种 算术 运算 ， 即 加 法 、 减 法 和 乘法 运算 ， 
做 除法 只 能 被 2 "形式 的 数 除 。 那 么 ， 这 样 的 计算 机 用 什么 
方法 做 任意 数 的 除法 呢 ? 

数 b 除 以 数 a 就 相当 于 求解 方程 4x 二 5。 因 为 计算 机 能 狼 
以 或 除 以 2" 形 式 的 数 ， 所 以 ， 我 们 假定 1/2<a<1 ( 杰 
则 在 方程 a*=5 两 边 同 乘 以 或 同 除 以 2 的 适当 每 指数 ) ， 把 
”方程 写成 

x*=( 1—a)zx+b (10) 
的 形式 。 取 x, ==5 作 为 方程 的 第 一 次 近似 解 ， 解 的 误差 用 a 
素 未 ， 也 就 是 说 ， 假 设 2, 十 ci 一 5/a， 则 从 方程 (10 ) 我 们 
得 到 
Xi 十 ai 一 ( 工 一 a) (2 十 al ) 二 5=( 工 一 z)zi 

+b+( 1-—a)a (11) 


思 为 172 <a< 1， 所 以 0< 1 一 “< 去 。 
由 于 因子 工 一 < 相当 地 小 ， 因 此 ， 略 掉 方程 (11 ) 右边 不 超 
过 去; 的 ( 1 一 a ) ec 项， 我 们 得 到 


3%1 十 aiS( 工 一 0)2: +b 
数 *s 二 (1 一 a ) x1 十 4 将 作为 * 的 下 一 次 近似 值 。 
再 假定 x, 的 误差 用 us 表示 ， 即 是 *, 十 as 二 b/a， 则 从 
方程 ( 10 ) 我 们 又 得 到 
xz 十 os 一 ( 工 一 a)x: 十 5 十 ( 工 一 a) ca: 
路 去 方程 右边 的 ( 1 一 a ) c: 项 ， 我们 又 得 到 一 个 近似 
等 式 
xs+tasT( 1—a)x,+b 
因此 ， 我 们 可 选取 xs 二 《1 一 a ) x*, 十 5 作为 下 一 次 近似 值 。 
用 同样 方法 ， 我 们 可 得 到 再 下 一 次 近似 值 
4 二 《1 一 a ) *s 十 5 等 等 。 一 般 地 ， 从 公式 
Xnti= (1—a)*nt+b (12) 
中 逐次 计算 x*1，*，，…… ，*r， 最 后 就 融 近 5/a 的 准确 什 。 
然而 ， 在 这 个 公式 里 ， 仅 仅 用 了 加 法 、 减 法 和 乘法 运算 ， 
这 就 是 说 ， 计 算 机 可 通过 这 个 公式 计算 5/a 的 值 。 
上 面 描述 的 做 除法 的 方法 ， 实 际 上 是 根据 无 穷 项 递减 等 
比 级 数 求 和 公式 。 实 际 上 ， 分 数 b/a 可 写成 


和 一 一生， ， 结 合 前 面 的 公式 ， 我 们 得 到 


T= yt (17e) 二 EC 工 一) 二 


十 5( 工 一 0 ) "+ (15) 
假定 级 数 的 前 4 项 和 为 x,，， 即 : 
18. 


xn 一 六 十 5( 1 —a) t+b( 1—a)"-! 
明显 地 
Xi 一 5 十 5 (1 工 一 a) + +h (1 —a)" tb (1—a)™ 
一 5 十 (1 一 ao)[b+D(1I 一 a) 十 … 
二 (1 一 oa)"1] 一 5 十 (1 一 oa)xn 
这 个 公式 与 公式 ( 12) 恰 好 相同 ， 因 此 ， 我 们 可 用 公式 (13) 
中 前 ?项 和 式 的 值 作 为 分 数 b ea 的 近似 值 **， 随 着 2 增 大 ， 


这 个 和 式 就 珊 近 整个 级 数 的 和 。 因 为 寺 <a< 1 且 0<1 一 e 


所 当 ， 故 级 数 (13 ) 是 递 大 的 。 
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5 . 逐次 逼近 法 求 平方 根 


+ 现在， 我们 来 讨论 怎样 应 用 逐次 近 近 法 求 平方 根 。 在 中 
学 , :我 们 学 了 由 某 个 已 知 数 的 平方 根 求 另外 一 些 十 进 制 数 平 
方 根 的 方法 ， 实 际 上 ， 这 可 看 成 是 逐次 逼近 答案 的 方法 。 然 
而 ,这 种 方法 相当 复杂 ,学 生 们 常常 机 械 地 使 用 了 它 ， 也 还 不 
完全 了 解 它 的 计算 原理 。 这 里 ， 我 们 要 介绍 另 一 种 古代 巴 比 
伦 使 用 的 方法 ， 这 种 方法 ， 古 希腊 的 几何 学 家 希 洛 也 曾 用 
过 ， 后 来 ， 它 被 人 们 渐渐 忘记 了 。 但 是 ， 现 在 在 电子 计算 机 
上 还 时 常用 这 种 方法 求 平方 根 。 

例如 ,假设 我 们 要 求 28 的 平方 根 , 首 先 选择 平方 根 的 某 个 
近似 值 , 比如 取 z: 一 5， 假定 它 的 误差 用 ac: 表示 ， 即 28 一 
%: 十 a 1， 为 求 出 xc:， 对 等 式 两 边 平方 ， 得 ， 

28 一 25 十 10a ,十 ai: 
即 
.+1001— 3=0 (14) 

于 是 ， 我 们 得 到 一 个 关于 a 的 二 次 方程 ， 解 此 方程， 得 到 
0 一 一 5 土 28， 因 此 ， 要 求 c: 的 精确 解 就 必须 计算 v 弱 。 
这 似乎 成 了 一 种 自身 的 恶性 循环 ， 要求 V28 就 必须 计算 1， 
要 计算 a 又 必须 考虑 28。 

下 面 的 论述 给 我 们 解决 了 这 个 问题 .近似 值 *,=S 的 误 
副 w, 是 小 于 工 的 ， 当 然 wc: 更 小 ， 因 此 ， 在 求 c: 时 ， 我 们 咯 
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方程 (14 ) 中 的 较 小 项 a，*， 则 得 到 一 个 关于 a 1 的 近似 方 
程 10a1 一 3 0， 由 此 ai 多 0.3。 
这 样 ， 我 们 就 求 出 了 误差 项 ,的 近似 值 ， 由 于 28=5 
人 
xz 一 5 十 0.3 一 
二 
就 是 说 ， 我 们 假定 近似 值 x 一 5.3 的 误差 为 4,， 令 V36=Xa 
十 q:， 等 式 两 边 平方 并 略 去 最 小 项 ac: ,我 们 得 到 28 挟 zs8 十 
2x:a:。 因 此 ， 
28 一 23 2 
22: 
此 式 表明 ,w38 的 第 三 次 近似 值 是 
tt 
出 zx: 一 5.3 代 入 上 式 我 们 得 到 za 一 5.2915…… 。 同 样 ， 由 近 
似 值 zs 一 5.2915， 又 可 得 到 下 一 次 近似 值 
28 十 2:” _ 
2xs 
一 般 来 说 ， 如 果 我 们 已 经 求 出 了 w28 的 近似 值 z*， 则 下 一 次 
近似 值 是 ， 


Q2: 各 


2 
ee cs). 


因此 ， 我 们 每 计算 一 次 ， 痢 得 出 一 个 比 前 次 更 精确 的 28 的 : 
近似 值 . 当 x+ 1 与 *» 的 差 小 于 给 定 精 度 要 求 时 ,就 停止 计算 。 
例如 ， 我 们 计算 w28 的 值 ， 精 确 到 0.0001， 四 次 近似 值 就 够 
了 , 即 取 v28==5.2915( 实 际 上 ,Xs 二 5.2915,x4 二 5.2915 ) 。 

上 述 方法 ， 同 样 可 用 来 求 任何 正 数 的 平方 根 。 因 此 ， 妆 
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Xn+1 一 


我 们 计算 ve 的 值 时 ， 先 选择 某 个 初次 近似 值 >， 然后 殖 借 
助 于 公式 


xngl 一 2 (16) 
逐次 计算 下 一 次 近似 值 。 


公式 (16 ) 还 可 以 由 与 求 V28 不 同 的 办 法 得 到 。 假 设 我 
已 经 求 出 了 Va 的 第 f 次 近似 信 xs， 因 为 V3 二/ . 2 2 


由 此 可 见 ，wWa 是 数 x， 所 以， 我 们 
取 这 个 几何 平均 值 的 近似 值 作为 < 的 算术 平均 值 。 
即 取 


2 
zi 去 (2m 十 汉 ) = -2 十 xm” 


2xn 
这 正好 就 是 公式 ( 16 ) 。 
因此 ,上 述 求 平 方 根 的 逼近 方法 在 于 : 每 一 步 都 以 数 x。 


和 -的 算术 平均 值 来 代 葵 它 的 几何 平均 值 。 


现在 ， 我 们 讨论 ， 逐 次 逼近 方法 用 于 求 平方 根 是 否 一 定 
能 得 到 正确 答案 ， 也 就 是 说 ， 是 出 现 阿 基 里 斯 和 乌龟 比赛 的 
情形 呢 ? 还 是 出 现 阿 基 里 斯 史 在 券 羊 后 面 的 情形 ( 数学 家 们 
称 第 一 种 情形 为 收 和 但， 第 二 种 详 形 为 发 散 ) ? 我 们 将 证 明 ， 
求 平方 根 的 方 法 是 一 种 收 倒 药 方法 ， 总 是 能 得 到 所 求 的 
结果 。 

。 为 此 ， 我 们 来 计算 两 次 逐次 逼近 的 误 差 a* 一 va 一 za 和 

Qnvi=V a—Xn+1 根据 公式 《16 ) ，an#1 可 写 为 : 

Qngi=VA—Xn#1=V A”— 以 
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Xn xn ata 

和 -六 汪 
但 由 于 x。* 一 2*。wa+a= (xX, 一 Va) 一 aas 
因此 


Cn 
as1 一 一 一 (17) 
Sn 


我 们 只 考虑 wa 的 正 近似 值 zs 的 情况 。 因 此 ， 由 公式 (17 了 
我 们 得 出 ， 所 有 的 误差 &2，Qs，…… ，aun…'… 都 是 负数 。 
换 句 话说 ， 从 第 二 次 通 近 以 后 的 所 有 近似 值 都 是 过 镜 近 似 
值 "， 而 第 一 次 近似 值 可 以 是 过 剩 的 或 不 足 的 。 

由 (17 ) 式 很 容易 证 明 ， 每 计算 一 次 近似 值 z,， 误 差 的 
绝对 值 至 少 减 少 为 原来 的 1 /2 。 实 际 上 ， 等 式 (17 ) 可 以 
写 为 ， 


0 sn—a i 1 Va 
nf Dx "ow (3 - 座 )o 
所 以 
-|L_va | . ~ 
Qnzl 2 2 Qn (13) 
但 是 ， 因 为 *.>> 0， 由 此 得 出 
工 一 va 一 上 
2 2x" <2 


昂 一 方面 ， 如 前 所 述 ， 当 ? 之 2 时 ， 有 xz 之 vay 所 以 


。 这 说 明 ， 算 术 平 均值 总 是 大 于 几何 平均 信 。 
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工 -va >0 


2 2xnr 
这 就 导出 不 等 式 

工 -wa | 一 上 

上 益 <5 (19) 
比较 关系 式 (18 ) 和 ( 19 ) 便 知 : eas < 二 |e， 


这 样 就 证 明了 我 们 的 论述 ， 每 一 次 融 近 计算 ， 识 差 的 纺 
对 值 少 于 前 次 误差 的 绝对 值 的 一 半 。 也 就 是 ，. 第 二 次 逼近 凡 
后 ， 误 差 的 绝对 值 就 少 于 最 初 误差 的 1 4， 第 三 次 逼近 以 
后 ， 就 少 于 1 8 ， 如 此 等 等 。 明 显 地 ， 随 着 n 的 增 火 ， 误 
差 的 绝对 值 cs 一 \76 一 和 将 不 断 减 少 ， 直 至 趋 于 零 . 即 是 说 ， 
x 趋向 Wa。 

现在 ， 我 们 来 讨论 ， 初 始 值 x, 的 选择 对 逼近 过 程 有 僻 
影响 ? 首先 ， 要 注意 到 前 面 已 经 证 明了 ， 无 论 按 择 怎样 的 初 
始 值 * ,， 都 不 影响 最 终 的 结果 ， 随 着 7n->o0, 误 牙 序 列 ,a ss 
as，…an… 趋 向 于 零 . 因 此 ， 只 要 给 定 了 所 需 的 精度 要 求 ， 
对 于 任何 初始 近似 值 x,， 我 们 都 会 得 到 Wa 的 满足 精 度 至 求 
的 同一 个 值 。 即 使 选择 一 个 差 值 很 大 的 初始 近似 值 ， 最 居 ， 
我 们 还 能 得 到 正确 的 结果 。 实 际 上 ， 经 过 10 次 逼近 以 后 ， 议 
差 的 绝对 值 至 少 会 碱 少 至 原来 的 1 A1000 (2'"=1024 沪 
1000 ) ， 经 过 40 次 逼近 以 后 ， 误 差 的 绝对 值 至 少 减 少 到 10 刀 
亿 分 之 一 (10!* ) 。 因 此 ， 当 我 们 计算 的 信 时 ， 如 果 取 
初始 值 * ,二 10°， 则 误差 a 1 守 10*。 而 | Q40 ! <10-., 也 就 
是 说 ， 开 始 计算 的 误差 大 约 是 工 百 万 ，40 次 逼近 以 后 ， 误 闭 
绝对 值 将 小 于 百 万 分 之 一 。 

然而 ， 初 始 近似 值 的 选择 要 影响 逼近 时 间 ， 如 果 初 始 近 
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似 值 选择 得 不 恰当 ， 一 般 需 要 很 长 时 间 ， 两 次 近似 值 xu+r :和 
x, 之 差 才能 小 于 给 定 精度 要 求 .选择 一 个 恰当 的 初始 近似 值 ， 
可 大 大 提高 逼近 速度 。 一 般 情 没 下 。 人 们 常常 是 从 平方 根 表 
中 取 - 一 个 适当 的 值 作为 初始 近似 值 ， 再 用 公式 


tx 
%> D7, (20) 


求 得 一 个 比较 精确 的 值 。 
这 种 求 平方 根 的 方法 特别 有 效 。 因 为 ， 当 x 逼近 Wa 时 s 
-误差 减少 的 速度 增加 得 相当 快 ， 这 是 因为 在 推导 不 等 式 


时 ,~ 我 们 用 去 代替 了 公式 (18 ) 中 


Qn+ 1 


1| 
< 了 lo 


4 |。 但 是 ，zn 接 近 vZae 时 ， 1-¥e 非常 小 ， 因 


呈送 地 主人 小 得 多 。 


我 们 还 可 以 更 准确 地 说 明 这 一 点 ， 为 此 ， 将 近似 值 * 的 
绝对 误差 | a | = | va 一 xz。 | 和 相对 误差 B* 一 起 讨论 。 相 
对 误差 就 是 绝对 误差 .| a* | 与 Va 的 精确 值 之 比值 ， 用 公式 
表示 即 是 


下放 | 


由 公式 (17 ) 可 得 到 


此 ， 


Cnr Qn 


因为 za， 所 以 


| an | 二 二 
BMS 2 


于 是 ， 相 对 误差 B， 满 足 不 等 式 
2 


Bari<iB, (21) 


例如 ， 如 果 近 似 值 *» 的 相对 误差 是 0.01， 则 x，4+ ;的 相对 误差 
不 超过 0.00005，x+, 则 不 超过 0 .0000000013。 

我 们 看 到 ， 随 着 计算 次 数 增加 ， 近 似 值 的 精确 度 以 极 高 
的 速率 不 断 增进 。 还 可 证 明 当 非常 接近 于 va 的 时 候 ， 每 一 
次 逼近 ， 可 以 使 精确 的 有 效 数 字 位 数 加 售 。 

例 ， 计算 v238， 精 度 要 求 为 0.00001 

从 平方 根 表 中 查 得 238 二 15.43 令 x*,=15.43， 再 用 


公式 求 xz: 
x 一 .15.432 十 238_ 
30.86 


现在 估计 答案 的 精确 度 。 因 为 15.43 的 误差 不 超过 0.01， 即 ， 


15.42725… 


人 一 0.01， 所 以 Bis 205<0.001， 因 此 


B:<-000L 一 0.0000005。 这 就 是 说 ， 近 似 值 x， 


的 绝对 误差 | c: | 不 超过 15.43X0.0000005<0.00001， 或 
者 说 ，238=15.42725 的 所 有 7 位 数字 部 是 精确 的 。 

如 果 我 们 想得到 14 位 精确 数字 ， 就 要 由 刚才 得 到 的 结果 
进行 第 三 次 允 近 。 然 而 ， 需 要 这 样 高 的 精度 是 非常 罕见 的 。 

最 后 ， 谈 一 下 逐次 允 近 法 的 特性 。 采 用 通常 的 方法 求 平 
方 根 时 ， 任 何 一 步 计 算 发 生 错 误 时 ， 就 会 使 后 面 的 计算 无 
效 。 但 使 用 逐次 逼近 法 就 完全 不 同 了 。 假 如 我 们 计 算 第 z 次 
通 近 的 近似 值 z* 时 发 生 错 误 ， 而 得 到 一 个 错误 结果 3*， 这 
时 ， 所 有 后 面 的 计算 就 相当 于 以 y, 作 为 计算 wa 的 初 始 近 似 
值 重 新 计算 。 但 是 ， 前 面 我 们 已 经 看 到 ， 无 论 选择 怎样 的 初 
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始 值 ， 用 逐次 豆 近 法 都 能 得 到 满足 精度 要 求 的 结果 。 因 此 ， 
即使 发 生 错 误 ， 计 算 中 的 误差 最 后 也 将 趋 于 零 ， 它 仅仅 给 我 
们 额外 增加 了 几 次 逼近 过 程 。 

由 于 逐次 逼近 法 的 这 种 特性 ， 计 算 可 从 较 低 精度 开始 ， 
对 指定 的 精度 只 需 最 后 几 次 逼近 。 这 样 就 大 大 缩短 了 计算 需 
要 的 时 间 。 
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6 .逐次 通 近 法 求 正 整 数 次 方 根 


上 面 讨论 的 求 平方 根 的 方法 ， 同 样 也 可 用 于 求 任何 正 整 
数 次 方 根 。 为 此 ， 我 们 需 用 公式 ” 。 
(x+a) 一 2 十 KE-1Q + (22) 
其 中 ， 省 略 符 号 表示 的 项 包含 kc*，a: 等 等 。 
现在 ， 我 们 来 证 明 这 个 公式 。 由 中 学 数学 课程 知道 ， 
(xz 十 a ) :一 22 十 2xa 十 aa2 
(x+a ) 一 22 十 3X2a 十 3xaa2 十 Qs 
这 起 每 式 ， 我 们 可 以 重 写 为 下 面 的 形式 
(x+a) :=x:+2x0 + (23) 
(xz 十 Ga ) *=x*+3xt0+ (24) 
芷 是 ， 这 就 证 明了 公式 (22 ) 出 k= 2 和 k= 3 的 情形 。 对 
公式 (24 ) 的 两 边 癌 乘 以 (x+a ) 得 到 
(x+a) 一 (X3 十 3x20 十 ……)》(x 十 G) 
销 去 括号 ， 我 们 将 质 到 ， 其 中 有 一 项 *4 不 包含 ac， 有 两 项 
包含 的 一 次 释 为 3xsa 和 zsa， 其 余 项 都 含有 a 的 二 次 和 更 
高 次 的 堵 。 所 以 ， 我 们 可 以 写 为 : 
(x+@ ) “=x*+3x0+x t=x+Ax 0 
十 四 二 (25) 
(与 前 面 一 样 ， 省 路 符号 表示 的 项 包含 a*、a :等 等 ) 


。 这 个 公式 杠 锯 二 项 式 定理 得 出 ， 但 我 们 不 要 求 汇 者 也 熟 严 人 交 个 定理 。. 
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这 样 就 证 明了 公式 (22 ) 中 k= 4 的 情形 ， 对 (25 ) 式 
施行 同样 的 方法 又 可 得 到 
(x#+a) 5 一 25 十 524Q 二 (26) 
明显 地 ， 我 们 同样 可 以 证 明 公式 〈 22 ) 对 任何 正 整 指数 
上 都 减 立 

现在 我 们 再 回 到 求 4 次 方 根 的 问题 上 来 ,这 里 & 是 任意 正 
整数 。 假 设 我 们 已 求 出 了 ka 的 某 一 近似 值 zx,，x: 的 误 问 
用 Qi 炎 示 ， 妈 是 说 ， 我 们 假设 有 Xi 二 ai 二 Ka, 那 人 么 (3 
十 ai ) “=a 由 (22) 式 ， 我 们 可 以 把 这 个 等 式 写 为 下 面 形式 ， 

xX1* 寺 kt io 二 =a 
其 中 ， 省 略 符号 走 示 的 项 包含 有 a,*，a1: 等 等 。 

如 果 选 择 的 近似 值 * ,非常 接近 Ka， 则 这 个 近似 值 的 
误差 1 是 非常 小 ， 因 此 ， 我 们 可 忽略 掉 含 有 a 高 次 等 的 那 
些 项 ， 得 到 下 面 的 近似 等 式 

1*+kt1i* ia Na 


册 此 得 出 ci 二 和 ar- ， 利 用 这 个 结果 ， 我 们 可 取 和 5 的 
下 次 近似 值 汶 : 


A md EA 


kx ZU 
同 理 ， 由 近似 什 x,， 又 可 求 出 下 一 次 近似 值 : 
i 4 十 (k 一 1 ) 2 
kes*™! 


一 般 地 ， 如 果 我 们 已 经 求 出 了 区 a 的 近似 值 z»， 则 下 一 次 
近似 值 可 由 公式 
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zw 一 .2 十 (大 一 1 ) 2s 
es 
求 得 。 

与 求 平方 根 的 方法 一 样 ,对 任何 初始 值 z, (要求 为 正 数 )， 
求解 方法 都 是 收敛 的 .换言之 ， 选 择 任何 x*,， 数 列 x1，*x，， 
Ds, Xn，…*… 趋 向 于 太 a。 通 近 过 程 要 计算 到 xn 和 xnzy 
都 满足 精度 要 求 时 为 正 。 

例 求 %970 的 值 ， 精 确 度 为 0.001 

在 近似 公式 (27 ) 中 ， 假 定 上 = 3 ， 即 


3 
GE 十 2Xw (28) 


这 里 ae 一 970， 令 x+ 一 10， 由 公式 (28 ) : 


970 十 2 xX103 ”2970 
eh ge en es 
3 X107 500 一 9.900 


970+ 2 x9.9s ”2910.60 
YX 一 Ey 一 
3X9.9: 294.03 一 一 9.899 


我 们 可 以 看 出 ， xz: 和 xs 都 在 精度 要 求 以 内 。 所 以 ， 以 970 精 
度 为 0.001 的 值 是 W970 二 9.899。 


7 .和 迭代 法 


前 面 讨论 的 所 有 例子 ， 都 是 求解 方程 一 般 方 法 的 具体 骞 
例 ， 这 种 方法 叫 和 迭代 法 或 逐次 逼近 法 ， 其 实质 如 下 

为 求解 方程 ( x ) = 0 ， 把 它 写 为 

z= (x) (29) 

的 形式 ， 然 后 选择 一 个 初始 近似 值 *:， 代 入 (29 ) 式 右 端 计 
算 ， 计 算 所 得 的 结果 *, 二 pq (*，) 作为 方程 根 的 第 二 次 近似 
值 。 一 般 地 ， 如 果 已 经 求 出 了 近似 值 x， 则 下 一 次 近似 亿 
Xn+1 可 由 公式 Xn+1 二 q(x ) 求 得 。 3 

假如 经 过 几 次 逼近 以 后 ，** 多 xz 里: 已 在 精度 要 求 之 肉 ， 
由 于 x+: 一 9(x), 这 就 是 说 xs (xn ) 也 满足 精度 要 求 ， 
即 x。 是 方程 x 二 pq《(%* ) 根 的 近似 值 。 

例如 ， 求 解 阿 基 里 斯 和 乌龟 的 问题 ， 我 们 把 方程 10* 一 4# 
一 1000 写 为 一 100 十 区 并 由 公式 zws: 一 100 十 名 
求 方程 根 的 近似 值 。 在 计算 机 做 除法 的 问题 中 ， 我 们 把 方 稳 
ay 一 5 写 为 Y 一 (1 一 a ) x 十 5， 并 由 公式 xn#1 二 (1 一 a)x。 十 
6 求 它 的 近似 值 。 最 后 ， 当 我 们 求 k 次 方 根 时 ， 把 方 程 x** ==@ 
化 成 


at(k—1)x* 


mm ti 


然后 ， 再 用 公式 
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x = + (k— 1 ) 2 
NF kn 
家 它 的 近似 值 。 这 里 再 给 出 一 个 用 迭代 法 解 比 较 复 杂 的 方程 
的 例子 。 
例 解 方程 10x 一 1 一 cos*= 二 0 (30) 
将 方程 ( 30 ) 写成 
_ 了 十 cosx 
C1) 
静 形 式 ， 选 择 一 个 初始 近似 全 。 比 如 取 x1 二 0， 代 入 方程 
( 31 ) 右 端 得 到 


1 十 cos0 
10 


我 们 用 它 作 为 所 求 根 的 第 二 次 近似 值 ， 并 将 *: 代 入 〈 31 ) 式 
右 端 ， 我 们 又 得 到 第 三 次 近似 值 


X2 一 一 0.2 


=_l1+cos0.2 ~ 1+0.98_ 
YXs 一 10 10 0.198 
然后 ， 再 求 出 
=_1+cos0.198 、 
0 0.198 


我 们 看 到 等 式 xs 一 2 已 满足 精度 0.001, 因为 24 一 


ee 这 就 是 说 ， 取 xs 一 0.198 是 方程 + 一 ces 


的 根 ， 具 育 0.001 精 这 。 
关于 选 代 法 解 方程 ， 有 以 下 几 个 问题 : 
图 用 先 代 法 得 到 的 序列 x:，x:，…xYn，… 是 否 一 定 收 全 

于 某 一 数 5? 


@ 如 果 没 限 , i 站 x 一 存在 ,是否 一 定 是 方程 + 一 p(x) 
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的 解 ? 

@@ 怎 样 才能 使 *1,，x:，…，x*。a… 迅 速 地 到 近 方程 x 二 9(%》 
的 根 ? « 

第 二 个 问题 很 容易 回答 ， 假如 zy X29 1 Xny … 通 近 
于 E， 考 虑 等 式 xn+1 二 p(x*，) ， 随 着 * 增 大 ， 等 式 左 边 带 近 
&， 等 式 右 边 通 近 于 pg(& )“， 因 此 ， 我 们 得 到 极限 值 &8 一 9 
(& ) ， 即 是 方程 x 二 q(x ) 的 根 。 

对 第 一 个 问题 ， 答 案 是 否定 的 。 例 如 :考虑 方程 ?一 10z 
一 2 ， 如 果 取 x, 二 1 ， 我 们 得 到 xz: 一 8， X%s 一 108 一 2， 
… 随 闭 " 增 大 ，xz 也 增 大 ， 并 不 趋向 于 任何 极限 值 。 另 一 方 
面 ， 如 果 我 们 把 方程 写 为 x<=1l0g( * 十 2 ) 的 形式 ， 通 近 过 
程 则 是 收敛 的 ， 并 且 经 过 三 次 和 逼近 以 后 ， 我 们 就 得 到 近似 人 
xX 一 0.38。 

所 以 ， 我 们 应 该 把 第 一 个 问题 换 为 如 下 问题 

什么 形式 的 函数 p( * ) 才能 保证 数 列 zi。，Xas，…?Za"o， 
收敛 ? | 

在 讨论 这 个 问题 之 前 ， 我 们 先 讨论 选 代 法 的 几何 意义 。! 


。 我 名 没 定 p(x ) 是 注 续 证 数 
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8 . 和 迭代 法 的 几何 意义 


显然 ， 求 方程 x 二 pg《* ) 的 根 E， 正 好 就 是 求 曲线 3 一 申 
《x* ) 与 直线 y= 二 x 的 交点 M 在 横 坐 标 上 的 值 。 设 有 某 一 初始 
值 *,( 图 2 ) ， 那 么 ， 点 Mi 的 坐标 MM， (21， q(x*,) ) 将 依 
赖 于 曲线 4 一 (x ) 。 过 这 点 作 一 水 平 线 与 直线 ;二 x 相交 于 
点 YNi(9 (xD)、 (xi))， 用 zz 表示 (xD)， 则 点 Ni 的 从 
标 为 N， (Za，Xa )， 然后 再 过 点 Ni 作 一 垂 线 ， 与 册 线 y=g 
(x* ) 相交 于 点 Ms (*:，9 (xz:))， 重 复 这 种 步骤 , 我 们 得 
到 ,与 直线 3 一 xz 的 交点 坐标 N:(zsxs)( 这 里 zs 一 pg(z:))， 
与 曲线 ;二 q(x ) 的 交点 坐标 M。( xs、9(x: ) ) ,如 此 等 等 。 
如 果 逼 近 过 程 是 收敛 的 ， 则 点 列 Mi, Mi, Ms, … 将 到 
近 我 们 所 求 的 直线 y=* 与 曲线 y= 二 pg( * ) 的 交点 。 

因此 ， 逐 次 带 近 法 的 几何 意义 就 是 ， 沿 着 一 条 折线 ， 向 
游 我 们 所 区 望 的 直线 和 曲线 的 交点 允 近 运动 ， 而 这 条 折线 的 


顶点 依次 出 现在 直线 和 曲线 上 ， 且 折线 段 是 垂直 和 水 平 交替 


出 现 的 (四 2a)。 

如 果 走 线 和 曲线 的 位 置 像 图 2 a 那样 ， 则 这 条 折 线 看 起 
来 就 像 一 架 梯 子 。 如 果 直 线 和 曲线 的 位 置 像 图 25 那 样 ， 则 
这 条 折线 君 起 来 就 像 一 条 螺 线 。 

有 时 候 遥 近 过 程 也 可 能 得 不 出 正确 结果 ( 阿 基 里 斯 和 羚 
羊 的 问题 就 是 如 此 ) ， 这 是 因为 点 列 Mi，M:，…，M。， 
… 不 是 朝 着 M 点 而 是 背离 M 点 运动 ( 图 8 ) ， 用 图 表示 就 是 . 
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梯子 ( 或 螺 线 ) 的 步 幅 变 得 越 来 越 大 。 

图 2 和 图 3 的 不 同 点 是 ， 过 直线 y 二 x 与 曲线 j 二 9 (*》 
的 交点 M 作 一 条 与 x 轴 的 类 角 为 135” 的 直线 ，`` 条 直线 过 
同 直线 y 一 * 把 空间 分 成 四 个 象限 。 如 果 在 点 M 邻 近 的 曲线 位 
于 左 象限 和 右 象限 ， 且 选取 初始 近似 值 也 在 这 个 邻近 区 域 
内 ， 则 选 代 过 程 收敛 。 反 之 ， 如 果 曲 线 位 于 上 象限 和 下 象 
腿 ,? 则 迭代 过 程 发 散 。 
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3 


(b) 一 、 


图 3 


然而 ， 使 用 这 条 规则 ， 首 先 要 描绘 出 函数 3=9 (zx ) 的 
粗略 图 形 ， 这 总 是 不 方便 的 ， 所 以 ， 还 得 考虑 另外 的 收敛 检 
验 法 ， 只 用 解析 法 而 不 需 作 任何 几何 图 形 就 能 确定 迭代 过 程 
收敛 或 发 散 。 这 种 检验 法 将 在 第 10 节 中 讨论 。 但 是 ， 我 们 先 
要 熟悉 压缩 映射 的 概念 。 
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9 .压缩 映射 


考虑 区 问 [ ae， 5 】 上 的 洱 数 y 二 9 (*)。 对 区 间 上 任何 
一 点 xo， 在 # 轴 上 都 有 一 个 对 应 点 4o， 即 Jo。 二 9 (*。)。 为 
了 作出 这 一 点 ， 我 们 可 以 过 xz 轴 上 的 点 ze 作 一 垂 线 与 函数 一 
钻 (z) 的 图 象 相交 ,然后 再 过 交点 作 水 平 线 与 ? 轴 相交 (图 4)。 


因此 ， 罗 数 y==q (x ) 给 出 了 区 河 ( a,，b ) 到 5 畦 的 一 种 映 
射 。 对 应 于 区 间 〔 co，4 ] 上 的 点 ， 被 陕 射 到 4 办 上 的 所 有 点 
的 集合 称 为 区 闻 [ a，4 ) 的 象 。 例 如 ， 商 4 三 x** 下 区 间 
[ 2，5 ] 的 象 是 区 间 [【 4，25 ) ， 在 同一 个 映射 下 区 间 
【一 1，6 ] 的 象 是 区 闻 【5 0，36 ] 。 可 以 证 明 ， 如 果 函 
数 ; 一 9 (2 ) 在 区 间 [ e，”] 上 连 著 ， 则 这 个 区 间 的 象 也 是 
y 轴 上 的 区 间 ， 并 且 如 果 画 数 3 二 w(x ) 是 一 个 单 35 递 增 函 
数 ， 则 区 间 < a，& ] 的 象 是 区 间 CC (a),，q(8)])， 如 
时 是 一 个 单调 递减 函数 , 则 象 是 区 间 (9()) ,9(e))。( 图 5 》 


我 们 讨论 了 区 间 (a,b ) 到 y 轴 的 映射 。 进 而 还 可 讨论 

它 到 x 轴 的 映射 ， 为 此 ， 区 间 映 射 到 4 轴 以 后 ， 将 3 轴 按 显 时 

针 方向 旋转 90” ， 即 是 区 间 〔 e,b ] 上 的 点 先 喘 射 到 4 轴 上 ， 

然后 ，3 贺 上 的 点 再 了 映 射 到 x 加 上 ， 用 这 种 方法 ， 函 数 g(x ) 

给 出 了 一 种 区 间 ( a，4 ] 到 xz 轴 的 映射 ， 我 们 把 这 种 喘 射 表 

示 为 :x 一 g(x ) 。 如 果 函 数 9〈2z ) 连续 ， 则 我 们 可 得 到 
一 个 x 轴 上 的 区 间 。 

也 可 能 出 现 这 样 的 情况 ， 区 间 [ao、# ] 的 象 C al5545] 

是 区 间 ( sa，5b ) 的 一 部 分 ， 例 如 : 在 映射 ?3 一 wz 十 1 下 ， 

33 


区 间 [ 0，4 ] 映射 到 它 的 部 分 区 间 [ 1，3 ), 在 这 种 
情况 下 ， 我 们 就 说 ，9 (*) 把 区 间 [a, $b 映射 到 它 的 
一 个 子 区 间 。 如 果 q(x*) 把 区 间 [ a,6 ] 映射 到 子 区 间 (c,，.、 
5，] ， 则 对 于 任何 子 区 间 ( o,5 ] 都 将 映射 到 子 区 间 (a,， 
b，) ， 特 别 ， 子 区 间 { a1， b1 通过 (x ) 又 映射 到 它 的 
一 个 子 区 间 ( a;，b, ， 在 同样 的 映射 下 ， 子 区 间 (4a;， 
b, ] 又 映射 到 它 的 子 区 间 ( as，5s ] ， 等 等 ,结果 我们 得 到 
一 区 间 系 : 

[a 6), (Car, pi]，… [ar brn), 
其 中 ， 每 一 个 区 闻 都 是 它 前 一 个 区 间 ( an，b。) 的 子 区 间 » 
由 此 [5 cn+i，bn+i 2] 是 映射 p (xz ) 下 [an，bn] 的 象 。 


例如 映射 :一 1 一 ;将 寺 -使 区 间 ( 0，4 ] 映射 到 它 


的 子 区 间 [ 寺 ， 号 ] ， 对 区 间 ( 于， 人 实施 同样 的 映射 
我 们 得 到 区 间 [ 总 ， 扩 ) 等 等 .每 一 个 区 间 都 逐次 包含 在 它 


前 一 区 间 里 。 

可 能 有 两 种 情况 ， 或 者 对 所 有 区 间 [ as，b ] 都 有 一 
个 共同 的 区 间 〔 c，d ] ， 或 者 这 些 区 间 只 有 一 个 公共 点 8， 
对 于 后 者 ， 就 说 区 间 系 收缩 到 一 点 &。 

下 面 , 我 们 将 用 数学 式 子 来 描述 区 间 系 [a,b3)，, [ary 
bs ) ，… [an，ba)… 收 缩 到 一 点 的 条 件 ， 为 此 ， 我 们 引进 
压缩 这 一 重要 概念 。 如 果 区 间 [ a，b ) 上 任意 两 点 间 的 距离 


通过 映射 9 (* ) 至 少 减少 到 二 ( M> 1 ) ， 则 9 (xz ) 是 把 区 


间 C4，5 ] 映射 到 它 子 区 间 〔 ak，5 ) 的 一 个 压缩 映射 。 因 
为 *1 和 x, 两 点 间 的 距离 为 1x:—*1 |， 所 以 压缩 映 射 条 件 可 
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由 公式 表示 如 下 ， 
如 果 存 在 数 4( 0 <4< 1 ) ， 使 得 对 区 间 〔o，8 ] 中 任 
意 两 点 1,，x，， 不 等 式 


lo (x —@ (x.) | <alx,—z, | (32) 
成 立 (4= 瑟 ) ， 则 称 映射 8 (x ) 是 区 间 [ o，5 ] 上 的 一 
个 压缩 映射 。 


区 间 Ca，56 ] 的 任意 子 区 间 ( cy ] 的 长 度 通 过 压缩 映 
射 9 (+ ) 至 少 减 少 到 M= 荆 。 实 际 上 ， 令 cl，d ] 是 区 间 


(Cc，4) 的 象 ， 则 c: 和 4 是 区 间 [c，d] 上 某 些 点 2 和 zs 
的 象 : 
ci 一 9 (2 ) di=9(zs) 
然而 
Ia-eol=}9(*) ~9 (x ) | <a |x,—x, | 
因为 点 *, 和 x, 在 区 间 ( c，4 ) 上 , 则 它们 之 间 的 距离 | *, 一 
%! | 小 于 区 间 Cc，4d ) 的 长 度 | d 一 c | ， 所 以 
lad—e | <g | d 一 c| 
这 就 证 明了 我 们 的 论述 。 
现在 ， 我 们 可 以 用 数学 式 子 来 描述 对 区 间 [ a,6 ] 逐次 
使 用 映射 p(* ) 得 到 的 区 间 系 Cal,，b1)，[ ca， 6517), 
【4an，bn】… 收 缩 到 一 点 的 条 件 ， 那 就 是 : 
如 果 [a，4& ) 到 它 子 区 间 5 oi bi ] 的 映射 p(x) 是 
一 个 压缩 映射 ， 则 区 间 系 〔 a,， bi ), Cas,s bs, 
[Can bn ] » … 将 收缩 到 区 间 C a, b) 上 的 一 点 &。 
实际 上 ， 由 于 9 (> ) 是 压缩 映射 对 任何 x 都 有 ， 
| bm 一 an | <al bm- 一 nl ! 


疝 还 ， 
1 05- 一 on- | <a 15-: 一 cv-z | 
因此 
| bn—an | <a’ | bs-2—an-2 | 
以 此 类 推 ， 我 们 得 到 
12 一 an | Sa | 58—al 
因为 0 之 1< 1， 所 以 ， 当 w 趋 于 无 穷 时 , 数 列 9，4*，…4"， 
… 趋 于 零 ， 并 且 区 间 Can，6b，) 的 长 度 | pb 一 as | 也 扑 于 
零 , 因 此 ， 不 存在 属于 所 有 区 间 (as，b，) 的 一 个 子 区 间 。 
所 以 ， 区 间 系 | 
Ca, by), Ca pb]，…[a bn) 
履 第 到 一 点 。 
设 后 ， 我 们 讨论 映射 9(% ) 。 不 等 式 (32): | 9(z: ) 
一 9《%1) | <4g | x 一 *1 | 对 任意 一 对 数 x ,和 x 前 是 满足 
的 ， 这 样 的 映射 是 整个 数 轴 上 总 一 种 压缩 ， 在 这 样 的 喘 射 
9( :2 下， 我 们 证 明 存 在 一 个 压缩 区 间 。 因 为 (32 ) 式 对 任 
意 两 点 1 和 x 都 满足 ， 所 以 只 要 证 明 存 在 一 个 区 间 被 p(x ) 
映射 到 它 自己 的 区 间 里 就 行 了 。 任 取 一 数 a, 合 b=qp (a )， 
再 选取 4 <1， 并 使 6 之 4q1， 令 


[2 一 a | 

工 一 4: 
我 们 将 汪 明 ， 区 闻 [4 一 了 ，a 十 了 ] 被 8 (2 ) 映射 到 它 的 一 
个 子 区 间 里 。 实 路 上 ， 令 ?是 区 闻 上 任 一 点 ， 则 |] x 一 a | < 


R， 由 于 不 等 式 (32 ) 得 ; 
lo (x)—5|=|9(x*)—p(a) | 
<a 12 一 4| <aR 

因而 
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19(xz) 一 ca =19(2z) 一 ?十 0 一 0 | 委 

19(x) 一 上 | 十 12 一 a 上 委 

1 一 a| 十 4R 一 4R 十 ( 工 一 4 ) 

一 (1 十 4 一 44 ) R< 
这 就 证 明了 区 间 [ a 一 R，a 十 R ) 上 的 任意 点 就 是 被 9 ( *) 
映射 成 同一 区 间 内 的 点 。 所 以 ， 映 射 p(* ) 压缩 了 区 间 (a 
—R, a+R)。 
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10. 迭代 法 与 压缩 映射 


现在 ， 我 们 再 回 到 迭代 法 上 来 ， 选 代 法 常常 用 于 求解 
* 二 q(x*) 类 型 的 方程 ， 如 果 E 是 方程 的 一 个 根 ， 则 & 二 q(5)， 
且 喘 射 *>q(x) 使 E 点 固定 在 它 所 在 的 地 方 。 因 此 ， 求 解 方 
程 * 二 p(x) 的 问题 就 相当 于 求 映射 p(*) 不 动 点 的 问题 。 

如 果 了 映射 p(*) 在 区 间 [a，&82 上 是 一 种 压缩 , 则 在 此 区 
间 上 总 有 一 个 不 动 点 。 为 求 出 这 一 点 ， 我 们 对 区 间 (c* 纪 逐 
次 使 用 映射 pg(x) 得 到 一 组 区 间 系 

[a1, b1) > [| 

因为 9(*) 是 区 间 (a，4b) 上 的 压缩 映射 ， 故 存在 唯一 点 8 属 
于 所 有 区 间 [as，bs】， 这 点 就 是 p(*) 的 一 个 不 动 点 。 

实际 上 ， 映 射 p(x*) 把 任意 区 间 [an，6") 映射 到 子 区 
伺 [antt， ba+1] 。 所 以 ， 区 间 〔cn， bn) 上 任意 一 点 x 的 
象 p(x) 都 在 子 区 间 (ant1，bar1] 上 ， 当 然 ， 也 在 区 间 
Lan，6n] 内。 因为 点 6 属于 所 有 区 间 [an，bn】 ， 它 的 象 
了 (5) 必然 也 属于 所 有 这 些 区 间 。 但 属于 所 有 区 间 的 仅仅 只 
有 一 点 5， 所 以 ，q(E) 二 E。 即 是 说 8 是 映射 q(x) 的 一 个 不 动 
点 。 

”这 样 ， 对 于 区 间 [a，5) 上 的 压缩 映射 ， 在 此 区 间 上 总 
有 一 个 不 动 点 ， 且 这 点 是 唯一 的 。 实 际 上 ， 若 存在 另外 的 不 
动 点 mw， m 一 9(9)， 则 有 不 等 式 

1m 一 5 一 19(m) 一 (5) | <a 1m 一 8 
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因为 0<3 鼠 1， 所以， 仅 当 | n 一 5&1 = 0 时 ,不 等 式 满足 r 
即 1 一 E。 

现在 ， 我 们 系统 地 阐述 移 代 法 收 全 的 必要 条 件 。 

假设 函数 q(x) 是 区 间 fa，5】 上 的 压缩 映射 。 则 对 区 
间 Ca, 5】 上 的 任意 点 %2， 由 Xa+1 二 (Xn) 得 到 的 数列 *。， 
X19%s，…， Xn， ""… 收 伍 于 这 个 区 间 上 方程 x 二 q(x*) 的 根 &。 

实际 上 ， 令 (an,，5.】] (n= 二 1，2，… ) 是 在 [a,6)】 
上 依次 应 用 映射 g(x) 得 到 的 一 系列 区 间 。 因 为 点 ze 在 区 间 : 
[a，5】 上 ， 它 的 象 x1 二 (x6) 在 区 间 [a1，81】 上 ， 区 间 
[as，651】 上 的 点 x 的 象 x, 二 q(%1) 在 区 间 [4s，42) 上 ， 
如 此 等 等 。 于 是 ， 对 任意 的 x?， 点 *, 都 在 区 间 [an，6n】 上。 
由 于 区 间 [an，bm] 的 长 度 随 着 n 增 大 而 趋 于 零 。 所 以 ， 点 
列 4,，xz:，…，Yn，… 逼 近 于 这 个 区 间 上 的 公共 点 8o 

上 面 论述 表明 ， 区 间 〔o，2] 上 的 任意 点 xz， 都 可 作为 
初始 点 。 

现在 ， 我 们 再 来 求 点 列 Y，Yik，…，Xxa，… 有 逼近 & 的 速 
度 。 因 为 8 一 9(8)， 所 以 对 区 间 [e，5] 上 的 任意 点 c， 过 - 

19(c) 一 5 = |9()—o(E)|<alc-el (33) 
把 不 等 式 ( 35 ) 应 用 于 点 列 Ye，2i，Y:，…，yYn …， 因 六 
sn 二 pq(%»-1)， 由 此 得 出 

[xr—E |= |9(xr-1)~—E | <q |xr-1—E | 
然而 对 任何 x 有 

1xa—E | <a xi 一 | <g’ 12- 一 | < 

<a" | xzo 一 | 

因此 ， 确 车 x 增 大 ， 误 差 | x* 一 上 | 减少 的 速度 至 少 和 公 比 为 
4 的 几何 级 数 碱 少 的 速度 一 样 快 。 

我 们 举 一 些 例子 说 明 怎样 应 用 上 面 证 明 的 条 件 。 
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例 1 迭代 法 能 用 于 求解 方程 :二 于 zs 吧 ? 《34) 


y A+ 
此 例 中 ，@9(x) 一 tHe 对 任意 z, 和 x: 有 


lot) -9 1 = 一 二 过 
EX x 二 | 
OIC 
i ls:+x*2| 
加 | 4 十 252 [|。| 4 十 22 | 
应 用 算术 平均 值 和 几何 平均 值 的 不 等 式 ， 我 们 得 到 


4 十 2: 
4 


* | *,—*x, | 


1z1= 二 /< 
Ixit+x2 | < 1x 1 十 |x: I< 全 4 十 2 


= X12 十 X%22 X%12 十 X%22 2X12X 8 
二 2 二 人 必 2 二 全 于 十 开放 


= 4 44.)( 4 +x:*) 


于 是 我 们 证 明了 对 任意 的 *, 和 x,， 不 等 式 
% 1 十 *> < 二 
(4 十 Yi12)(4 十 Yaz2)、8 


成 立 。 所 以 
| 9(xs)—p(2) [< 二 1z: 一 2 1 


这 就 是 说 ， 映 射 9(*) 是 整个 效 抽 上 的 一 种 压缩 。 
我 们 已 经 知道 ， 在 这 种 情况 下 ， 存 在 一 个 区 各 被 压缩 映 
射 到 它 自己 的 区 问 里 。 为 求 出 得 数 ， 我 们 取 e= 0， 则 映射 
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9(z) 把 点 o= 0 映射 成 点 2 一 下。 此 外 ， 在 这 里 4 一 二 ， 我 他 


 _ 1 |b-al_lim_l 1 
取 41 一 十, 并 且 R 为 二 8| 一 吉 , 则 (一 告 s 吉 ) 是 被 9(*) 
映射 到 它 自 己 的 区 间 。 所 以 ， 在 此 区 间 里 ， 存 在 一 个 不 动 
点 , 即 是 方程 (34) 的 根 。 为 求 出 此 根 , 任 取 区 间 〔 一 去， 二 


上 一 点 ze= 0 作为 初始 点 ， 使 用 选 代 法 得 到 ， 
1 


“= 证 =0.25 
en er a ee 

“1 T0251 4.0605™ 0°2461 
= Ee 

Xa TF0. 23461 4.0605™ 0*2463 


Pr We 
4+0.2463: 4.0605 


所 以 ， 在 精度 0.0001 内 ， 有 xs 一 z4。 由 此 可 见 ， 取 精度 为 
0.0001 时 ， 方 程 ( 34 ) 在 区 间 [一 各 二 上 的 根 是 0.2463。 
因为 映射 p(x) 是 整个 数 轴 上 一 种 压缩 ， 所 以 方程 (34) 没 
有 另外 的 根 。 

例 2 能 用 逐次 盈 近 法 求解 区 间 [一 1，8 ] 上 的 方程 
二 1 十 3/x 胸 ? , | 

这 里 ，qP(*) 二 1 十 2Yx。 因 为 9( 一 1)= 0，9(8)= 3， 
所 以 ，q(*) 沁 区 间 (一 1，8) 映射 到 它 自 己 的 区 间 里 。 
然而 ， 在 这 个 区 间 上 ，9(x) 不 是 一 种 压缩 。 例 如， 假使 取 
xi 一 一 0.008，x: 一 0.008， 则 有 


二 0.2463 
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| p(x2)—q(x1) | = {3/0.008—.3-0.088 | 
=0.4> jx,—zx: | 
在 例 1 中， 我们 证 明 映 射 q(x*) 是 一 种 压缩 用 了 不 等 式 
本 < 中 *。 现 在 ， 我 们 再 介绍 几 个 不 等 式 ， 这 些 不 等 式 党 


党 被 用 来 证 明 某 种 映射 是 压缩 。 
证 明 对 xz> 0， 不 等 式 
sinx<x<tgx (35) 
成 立 。 为 此 ， 注 意 到 扇形 04B 的 面积 Sos4 (图 6 ) 介 于 加 
心 角 为 * 的 三 角形 04C 和 0AT 的 面积 之 间 ， 即 
SAo4a<Sseero04B<SAor4 
但 由 于 


《 R 为 加 的 半径 ) , 且 扁 形 0 4AB 的 面积 是 宇 。- 
因此 > 
约 去 B ， 我 们 得 到 不 等 式 ( 35 ) 。 


从 不 等 式 ( 35 ) 还 可 得 出 : 
对 0<x< 1 有 x<arcsinx 2 
对 x 汪 0 有 Xx>arctgr 
另外 还 有 不 等 式 图 6 
er>1+x (x>0) 
In(1+x*)<x (0<x<1) 
要 证 明 它 们 稍微 困难 些 。 
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例 3 ”用 迭代 法 求 出 方程 
x=1 + arctgr (36) 
的 根 。 
因为 对 所 有 的 x 值 ， 有 1 十 于 aretgx> 0， 所 以 方程 多 
有 正 根 。 这 里 
9(z) 一 工 十 二 arctgx (37) 
因此 
19(x:) 一 (xD 一 | (1 + 于 aretgz:) 


一 (1 十 二 aretgxi) | = | arctgy: 一 arctgxi 上 


但 对 x: 之 0，x: 过 0 有 


arctgx2 —arctgx1 =arctg 生 2 
所 以 
| p(x*2)—q9(*1) =z 直 re A 
(2)—m(#1) | 生生 


了 | 和 2 一 %1 
< 站 生 < 吉 :一 * 


由 此 可 见 ， 喘 姻 qg(*) 是 举 轴 [ 0，<。 ) 上 一 种 压缩 ， 它 把 区 
间 【0，w 8] 了 欢 射 到 它 自己 的 子 区 间 〔【 1， 工 十 车 ] 上 。 


所 以 ， 在 区 间 〔 1， 1+$】 上 有 方程 ( 36 ) 唯一 的 根 。 尖 
求 此 根 ， 取 x1 二 1， 则 
X2 一 1 + 二 arctg1 一 1 + 名和 Y1.39 
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xs 一 1 + 入 arctg1.39=1.474 


2 一 1 十 二 arctg1.474=1.487 


xs=1 + 二 arctg1.487 一 1.489 
Xe 一 了 十 去 arctg1.489 一 1.490 


YXYr 一 工 + 十 寺 aretg1.490=1.490 


可 见 ，xe 一 #r 一 1.490 已 满足 0.001 的 精度 。 这 就 是 说 ， 方 程 
具有 这 个 精度 的 根 是 1.490。 因 为 映射 w(*) 是 整个 半 轴 
0 x< 吕 上 的 压缩 ， 所 以 方程 ( 36 ) 没有 另外 的 根 。 
如 果 遇 到 一 个 方程 z=qp(*) 不 能 用 迭代 法 求解 ， 可 以 通 
过 变换 ， 化 成 可 用 迭代 法 求解 的 方程 。 例 如 ， 我 们 考察 方程 
X 一 Xs 一 2 (38) 
这 里 ， 因 为 pg(1)= 一 1<1，9(2)=6>2。 所 以 在 区 间 
[1，2] 上， 这 方程 有 一 个 很 。 但 是 ， 映 射 * :一 2 在 这 个 
区 间 上 却 不 是 一 种 压缩 ， 因 为 它 不 把 区 间 ( 1，2 ) 映射 到 
它 自己 的 子 区 间 。 我 们 把 方程 ( 38 ) 重新 写 为 
Xx 二 Xx 十 2 
的 形式 ， 令 Hz) 一 &2 士 2， 我 们 有 
| Wx2)—w(21) | = | sx 十 2 一 x 二 2 | 
2 一 XL 
YCxst2) + ta 2) + Yt2)i 
在 区 间 【1，2] 上， 有 x*i 之 1， *s 之 1， 所 以 
1 


[yx2)—Y(z1) | <ays 


| 
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这 就 证 明了 映射 W(x) 是 区 间 [1，2] 上 的 一 种 压缩 。 取 
2 一 1， 应 用 迭代 法 得 : 

Xs=/ 3=1.442 

xs=8/3.442=1.510 

x4=8/3.510=1.520 

xs=8/3.520=1.521 

Xe 一 3/3.521 一 1.521 
因此 ， 在 精度 0.001 内 ，1.521 是 方程 ( 38 ) 在 区 间 (1，2) 
上 的 根 ， 方 程 没有 其 它 的 根 。 

我 们 罩 到 ， 对 原始 方程 作 些 适 当 变形 ， 就 可 化 成 能 够 用 : 
选 代 法 求解 的 形式 。 

这 里 讨论 的 和 迭代 法 收敛 性 检验 ， 使 用 不 是 很 方便 的 ， 因 
为 需要 证 明 很 复杂 的 不 等 式 。 下 面 ( 见 第 21 节 ) ， 我 们 将 还 
要 讨论 一 种 收敛 性 检验 定 吾 ， 用 它 证 明 选 代 法 的 收敛 性 就 非 
常 容 易 了 。 ， 


11. 藤 截 法 


迭代 法 是 求 方程 近似 解 最 基本 的 方法 之 一 ， 其 它 许多 求 
方程 近似 解 的 方法 都 是 迭代 法 的 特殊 情况 。 现 在 ， 我 们 讨论 
其 中 一 种 方法 ， 叫 弦 截 法 ( 试 位 规则 ) ， 或 者 叫 弦 位 法 ( 比 
例 法 ) 。 

假设 我 们 要 解 方 程 1(x) = 0 ,这 就 相当 于 求 函 数 y 二 f(x) 
的 图 像 与 z 轴 的 交点 。 假 定 函数 f(x) 连 续 ， 且 在 a 点 和 5 点 ， 
f(x) 的 值 有 不 同 符号 ， 则 在 区 间 (a，b】 上， 至 少 有 一 点 使 
函数 为 零 。 换 句 话 说， 函数 y=f(*) 的 图 像 在 区 间 (a，6】 
上 至 少 有 一 点 与 * 轴 相交 。 一 般 情况 下 ， 可 能 有 几 个 这 样 的 
点 (图 7 ) 。 因 此 ， 如 果 函 数 # 一 f(z) 在 区 间 [a,，b) 上 是 
单调 函数 ， 且 在 区 间 的 两 个 端 
点 的 值 有 相 异 符号 ， 则 函数 图 ” 

像 与 + 轴 仅 有 一 个 交点 &。 为 求 2 

出 这 一 点 的 近似 值 ， 用 弦 MN 0 bx 
代替 曲线 y==f(x) 在 区 间 [a， 

5b) 上 的 一 段 弧 ， 并 求 出 弦 与 * 

轴 的 交点 T (图 8 ) 。 图 7 5 

为 此 ， 考 虑 相似 三 角形 MM 7 和 NN ,PY， 由 三 角形 相似 
得 出 . 
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图 8 


但 从 图 8 由 可 以 填 出 : MT=a,~as, ZN 一 0 一 ay MM 
二 ~) NN=1(4)。 其 中 ，4 ,表示 绽 MN 与 * 轴 交点 的 
横 你 标 。 所 凡 


4 一 5 ba 
~tf(a) ff(5) 


解 这 个 方程 ， 我 们 得 到 


地 可 写成 
rs b--a 
a b fb) Tay (39) 
或 
ef rp bey (40) 


ai 就 是 方程 if) 二 0 在 点 a 和 点 b 之 闻 根 的 近似 信 。 

因为 f(a) 和 f(5) 的 符号 相 异 ， 所 以 有 两 种 可 能 ， 或 者 
f(a) 与 1(a,) 和 的 符号 相 异 ， 或 害 f(5) 与 (a4) 的 符号 相 异 。 如 
果 函 数 f(z) 在 点 4 和 ai 值 的 符号 相 异 ， 则 把 公式 (39 ) 应 用 


4 


于 区 间 (a，ai 】 上 ， 从 而 得 到 所 求 根 的 近似 值 


ai 一 4 
FeOD) 二 7 (41) 

如 果 是 另 一 种 情况 , 即 函 数 f(xz) 在 点 al 和 的 值 有 相 异 符号 ， 
册 把 公式 (40 ) 应 用 于 区 间 [ai，5]， 得 到 所 求 根 的 近 . 
似 值 


az 一 0 一 ai) 


ao 一 fa 和 (42) 


求 出 a 的 值 以 后 ， 继 续 把 公式 ( 39 ) 应 用 于 区 间 (a,as7 
上 (或 者 把 公式 (40) 应 用 于 区 间 (o: ,上 ) ， 就 能 求 出 下 一 
次 近似 值 as。 一 般 地 ， 如 果 近 似 值 x" 已 经 求 得 ， 下 一 次 近似 
信 可 用 公式 


二 Cn 一 G 
Gari Tan f(an) pra Fey (43) 
残害 
_ b—an 
| @n+i 一 0n fa) 了 六 一 TD (44) 
求 得 、 


我 们 得 出 了 ‘4: ) 和 (44 ) 两 个 公式 ， 使 用 时 应 该 用 哪 
个 公 臣 呢 ? 假设 曲线 是 向 上 四 的 ， 在 这 种 情况 下 ， 曲 线 上 的 
各 点 簿 这 与 函数 为 正 值 的 任 一 端 连接 ， 无 论 它 是 M 端 点 或 N 
端点 、 到 之 ， 和 如 果 曲 线 是 向 下 四 的 ， 则 点 应 该 与 函数 为 负 值 
的 那 一 端 连接 。 图 9 中 说 明了 函数 图 像 可 能 出 现 的 不 同位 
置 。 这 些 图 像 ， 给 出 了 下 面 论述 的 直观 意义 。 

假设 函数 fx) 在 区 各 [ a，86 ] 上 单调 连续 ， 函 数 的 止 止 
帖 在 区 同 . 上 不 改变 方向 ， 且 函数 在 区 间 两 个 端点 的 值 有 祖 异 
将 号 ， 则 选 定 正确 的 台 近 公式 ， 由 弦 截 法 得 出 的 点 列 收敛 于 - 
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方程 i(x) = 0 的 根 。 
反之 ， 如 果 选 择 的 公式 不 正确 ， 由 弦 截 法 得 到 的 点 a;， 
就 可 能 超出 区 间 ( a,，b ] ， 图 10 就 说 明了 这 种 情况 。 
上 述 的 弦 戴 法 正好 就 是 迭代 法 的 一 种 特殊 情况 。 假 设 函 
数 f*) 在 x 二 a 点 不 为 零 ， 在 这 种 情况 下 ， 方 程 f(x)== 0 就 相 
当 于 方程 
49 


一 xz 一 f(z) -2 | (45) 


71(x)—i(a) 
实际 上 ， 如 果 i(&)= 二 0， 则 
运 ER -es 
EE—i(E) 0 fa (46) 


相反 地 ， 如 有 果 E 二 a， 县 ( 46 ) 式 成 立 ， 则 1(&) 二 0。 
然而 ,方程 (45 ) 是 x 二 q(x) 的 形式 。 其 中 


X 一 0 af(x)—zxf(a) 
f(x)—f(a) f(x)—f(a) 


取 x。 二 5， 应 用 迭代 法 ， 我 们 将 得 出 一 组 数列 


外 (X) 一 2 一 人 2X)》 


Qi Ary ee any 
这 组 数列 与 用 弦 截 法 公式 
n+1l 一 0n 一 fan) per fa) 
得 出 的 数列 一 样 。 
作为 例子 ， 我 们 用 弦 截 法 求解 方程 
x 二 3x 一 1 二 0 (47) 


这 里 ，1(x) 二 x: 十 3x 一 1。 因为 1(0)= 一 1，f(1)=3， 所 
以 方程 在 区 间 ( 0，1 ] 上 至 少 有 一 根 。 如 果 我 们 作出 函数 
9 二 ** 十 3x% 一 1 的 图 像 ， 可 以 看 出 ， 它 在 区 间 ( 0，1 1] 上 
的 图 像 是 向 上 四 的 。 所 以 、 我 们 使 用 公式 ( 39 ) 。 按 照 这 个 
公式 ， 方 程 第 一 次 近似 根 是 


2 关 
A 0 9 


为 求 出 第 二 次 近似 值 ， 我 们 用 公 


一 1 nd 025 
A ne 一 1 一 3X SF0.23 


=0.31 
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于 是 
xs 一 1 一 3Xx 工 二 0.31 0.319 


3 十 0.040 
1。 一 0.319_- 
x4=1—3x3F0 0 322 
de 1 一 0.322 _ 
#6— 1— 83X30-0006™ 322 


所 以 ， 方 程 在 区 间 (0，1) 上 具有 精度 0.001 的 根 是 0.322。 
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12. 改进 的 弦 截 法 


人 
根 的 误差 值 按 等 比 几何 级 数 减少 。 有 一 种 改进 的 弦 截 法 ， 可 
使 收敛 速度 大 大 提高 。 使 用 弦 截 法 时 ， 每 一 次 计算 近似 人 都 
用 了 区 间 [ o，5 的 两 个 端点 值 之 一 。 我 们 可 以 用 前 两 次 计 
算得 出 的 近似 信 来 代替 区 条 Co, 4 ] 的 两 个 端点 信 ， 因 为 它 
们 更 接近 方程 的 根 。 

使 用 前 两 次 近似 值 计算 的 公式 形式 如 下 ( 图 lla ) ， 


anti=on—f(an) pea fy (48) 
这 里 ，a! 由 公式 (39 ) 计算 ，c: 由 公式 (41) 或 (42) 计 
算 ， 取 决 于 f(a)，f(8) 和 f(a1) 的 符号 而 定 ， 如果 f(a)<< 0， 
f(0) >0， 则 feat)< 0 时 选择 公式 (42 ) ，f(ai)> 0 时 选 
择 公式 (41 ) 。 
如 果 由 公式 ( 48 ) 计算 出 的 4:， 结 果 是 落 在 区 间 (a,8) 
之 外 ， 计 算 下 一 次 近似 值 时 ， 就 把 举 近 as 的 那个 区 间 端 点 值 
用 来 代替 as( 图 118 ) 。 
改进 弦 截 法 的 收敛 性 比 通常 弦 截 法 的 收敛 性 ， 结 果 要 好 
得 多 。 妈 ， 如 果 E 是 方程 1(#)== 0 的 根 ， 则 有 
lan+1 —8l< cian —A' (49) 


其 中 ， := lt Sl.6l8 
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作为 例子 ， 人 
X3 十 3x 一 1 一 0 
与 通常 弦 截 法 一 样 ， 第 一 次 近似 什 a, =0.25, “ 且 a; =0.31。 
下 一 次 近似 值 就 由 公式 (48 ) 计算 : 
Fa fy 0.31+0.040x 二 
一 0.3223 

我 们 有 #(0.3223) 二 0.0004。 显 然 ，*= 二 0,3223 作 为 方程 

的 根 ， 精 确 到 0.0001。 


4 一 0: 一 f(a:) 


13， 多 项 页 式 的 向 商 


借助 于 迭代 法 求解 方程 f(x)= 0， 很 大 程度 上 取决 于 怎 
样 把 方程 化 成 + 二 q(x) 的 形式 。 许 多 情况 下 ， 最 好 采 用 牛顿 
根据 微 商 ( 即 导数 ) 概念 提出 的 方法 。 这 一 节 里 ， 我 们 要 讨 
论 多 项 式 的 微 商 并 对 微 商 作 些 解 释 ， 使 我 们 能 用 牛顿 方法 去 
求解 代数 方程 ， 即 下 面 形 式 的 方程 

aoXx 十 aiXt1 十 … 十 ar 一 0 (50) 
令 

F(X)=aox*+arx* 1! 二 .+ax 
f(x) 是 一 个 多 项 式 。 考 虑 多 项 式 !(x 十 a)， 即 式 子 

ao (x 十 a)* 十 a1 (x 十 Qa)*-! 十 … 十 ax (51) 

如 果 消 去 (51 ) 式 中 的 括号 ， 我 们 发 现 ， 有 些 项 完全 不 
含有 ac， 有 些 项 含有 w 的 一 次 窜 ， 有 些 项 含有 au 的 二 次 寡 ， 如 
此 等 等 我 们 把 含 c 同 次 需 的 项 台 并 在 一 起 , 则 多 项 式 fx 十 a) 
可 表示 成 下 面 的 形式 
f(x 十 Q) 二 fo (Xx) 十 f(x)a 十 f(x%)a? 十 … 十 fe (x)a* (52) 
( 因为 多 项 式 f(x) 的 最 高 次 数 是 k 次 ， 所 以 展开 式 (52 ) 中 a 
的 最 高 次 数 也 为 8 ) 明显 地 ，fo (x)，f1(x)，…，fr(zx) 也 是 
关于 x 的 多 项 式 

例 取 f(x)=2x? 一 3x? 十 6x 一 1， 则 

f(x+o)=2(x+a):—3(x+a) +6(x+a)— 1 

一 2(X3 十 3x2a 十 3xa2 十 2) 一 3(z2 十 2zQ 十 2) 
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十 6(z 二 a) 一 下 一 (2xz2 一 3x2 十 6x 一 直 ) 
十 (6x* 一 6x 十 6)a+(6x 一 3)a: 十 2a* 
所 以 此 例 中 
ff(z) 一 2xs 一 3xz2 十 6xz 一 1 
f1(x*)=6x*—6x+ 6 
f(x)=6x—3 
fs (x)=2 
我 们 看 到 ，f。 (x) 与 f(x) 一 致 ， 这 决 不 是 偶然 的 巧合 。 如 果 
我 们 在 方程 ( 52 ) 中 令 a= 0， 便 得 到 f(x) 二 fo (x)。 
现在 ， 我 们 再 看 下 一 项 f1(*)a，a 的 系数 即 多 项 式 f(x) 
就 称 为 多 项 式 i(%) 的 微 商 。 例 如 多 项 式 2x* 一 3%* 十 6x 一 1 的 
微 商 是 6x: 一 6x 十 6 。 一 个 多 项 式 的 微 商 通常 记 作 # (x)。 
因此 ， 多 项 式 f(x) 的 微 商 /' (x) 就 是 多 项 式 t(x 十 a) 展 开 


式 中 关于 a 乘 竹 项 的 系数 。 
引用 上 面 的 记号 ， 我 们 可 以 把 ( 52 ) 式 重 写 为 
f(x+a)=f(x) + (x) ot (53) 


其 中 ， 省 略 符号 表示 的 项 包含 a:，a*，……，ax。 例 如 
2(x 十 ao) 一 3(X 十 a) 2 十 6(X 十 a) 一 工 
二 2x3 一 3%2 十 6x% 一 1 十 (6x* 一 6x 十 6 )a 十 … 
我 们 介绍 了 多 项 式 f(*) 微 商 的 概念 ， 现 在 再 讨论 怎样 计 
算 微 商 。 为 此 考虑 多 项 式 
f(x+a) 一 ao (x 二 a)* 二 ay (x 十 Qa)* 十 … 十 ar- (x 十 a)》 
十 ak 
对 其 中 每 一 项 者 用 展开 式 (x 十 q)"=x" 十 mx"” ig 十 … 代 换 
《 见 第 6 节 )， 我 们 得 到 
f(x+a)=ao (x*+kx*- lat.)+a tx*-! 
十 (一 1)Xt2a 十 …] 十 … 十 az-1(X 十 G) 十 ax 
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一 Go3Xx 十 CiX4 1 二 rata(taox*-! 
十 (一 1)ak3xt2 十 … 十 ac- 十 …… 
把 这 个 等 式 与 ( 53 ) 式 
f(x 二 Qa) 一 AZ) 十 af (2) 十 … 
比较 ， 我 们 得 出 ; 
多 项 式 
fz) 一 0oXt 十 aiXt 1 十 … 十 ak-iZ 十 ax (54) 
的 微 商 是 
f(x)=iaort-!1+(k—1)axt ?+ +a, (55) 
例如 ， 多 项 式 
f(xz) 一 637 十 8&Xs —3x:— 1 
的 微 商 是 
f' (x)=42x° +24x2—6x。 / 
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14. 牛顿 法 求 代数 方程 近似 解 


现在 ， 我 们 再 回 到 用 牛顿 法 ( 切线 法 ) 求 代数 方程 近似 
解 上 来 。 假 定 有 方程 
aoxx 二 aixt1 十 … 十 ax 一 0 (56) 
假设 已 经 以 某 种 方法 求 出 了 方程 根 的 一 个 近似 值 zx:， 我 们 将 
说 明 如 何 求 得 这 个 根 更 精确 的 值 。 设 近似 值 ?: 的 误差 是 aiy 
即 x1 十 a, 是 方程 (56 ) 的 根 ， 则 必 有 
ao (Nita) +a (tT) -1+'a:= 0 (57) 
也 就 是 说 
f(x1+ai)= 0 
这 里 1(*) 表 示 多 项 式 aox* 十 a1x*"! 十 … 十 as 
但 是 ， 按 照 公式 (31 ) 有 
jx1 二 1) =F(X1) of (#1) + 
其 中 点 表示 的 项 包含 ga1?，@1 5， aixs。 因此 ,为 确定 ciy 
我 们 有 方程 ‘ 
> fx 十 ai) 一 (xi) 十 aif (Xi)+t+ =0. (58) 
如 果 初 始 近 似 值 x 足够 精确 ， 则 误差 1 很 小 ， 这 样 ， 方 
程 ( 58 ) 中 省 略 的 项 相对 于 a, 是 非常 小 量 ， 略 去 这 些 项 ， 我 
全 得 到 一 个 决定 a 的 近似 方程 < 
f(x1) +af’ (x TO (59) 
由 此 得 出 
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~ f(x1) 
en) (60) 


所 以 方程 根 的 改进 值 就 是 


x = (61) 


然后 ， 我 们 再 改进 得 到 的 近似 值 ， 求 出 第 三 次 近似 值 


fxa) 
f(x,) 


一 般 地 ， 如 果 我 们 已 经 求 出 了 方程 根 的 第 n 次 近似 值 ， 则 下 
一 次 近似 值 计算 公式 是 
x Tx) 
YXn+1 一 Xn Fr Cxn) (62) 
公式 详细 形式 如 下 : 


Xntl Xn 


Xs=Xs 


GoXn* 十 @1Xn* 1! 十,… 十 Qs-1xn 咎 a* 
kaoxn*-!+(k—1)alr* 2 十 … 十 ax 
(63 ) 

如 果 近 似 值 z* 和 x+; 满 足 精 度 要 求 ，( 在 精度 要 求 范 转 
内 ) 也 就 求 出 了 根 的 近似 值 。 

上 面 讨论 的 求解 方法 应 归于 著 名 的 英国 数学 家 牛顿 
(NEWTON ) 。 

牛顿 方法 与 迭代 法 有 着 非常 紧密 的 联系 。 尤 其 是 ， 如 果 
函数 一 所 z) 和 3 一 (2) 没 有 共同 的 根 ， 则 方程 Kxz)= 0 就 
相当 于 方程 


x 荆 % 二 站 2 (64) 


对 这 个 方程 应 用 迭代 法 ， 我 们 得 出 一 列 数 x4，z4，…，Xa， 
…。 计 算 这 列 数 的 关系 式 
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zt (65) 
， 与 牛顿 法 公式 一 样 。 也 就 是 说 ， 牛 顿 法 在 于 把 方程 {x*)= 0 
写成 (64 ) 的 形式 ， 然 后 再 对 它 应 用 迭代 法 。 

例 用 和 牛顿 法 解 方程 x* 一 3x 一 5 = 0 要 求 精度 为 0.001， 
选 初始 近似 值 为 *, = 二 3 

因为 多 项 式 f(x)=x: 一 3x 一 5 的 微 商 是 1'(x*)=3x*:* 一 .3 
应 用 公式 ( 62 ) : 


| tnt tn 
所 以 
En 27 二 9 一 5 13_ 
人 27 一 了 24™2.46 
二 _14.89 一 7.38 一 5_ 二 
Xs 一 2.46 ~ ig.16—3 2.46 一 0.165=2.295 
一 2 295 一 12.088 一 6.885 一 5 一 一 = 
X4 一 2.295 15.801 二 3 2.295 一 0.016 一 2.279 
xu 一 2.279 一 也 .837 一 6.837 一 5 2? .279 


15.582 一 3 
我 们 看 到 ， 取 z4 一 zs 已 精确 到 0.001。 所 以 ,方程 2 一 3X 
一 5=0 精 确 到 0.001 的 根 等 于 2.279。 
在 第 6 节 中 介绍 的 近似 计算 整 次 方 根 的 方法 也 是 牛顿 方 
法 的 一 种 特殊 情况 。 实 际 上 ， 求 Ba 的 值 就 是 解 方程 
x*—a=0 
因为 多 项 式 *x* 一 a 的 微 商 是 tx*-!， 所 以 对 于 方程 
xx 一 4 一 0 
公式 (62 ) 为 ， 
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An* a_at(k—1)x,* 
了 汪汪 kxn*™! 
这 正好 就 是 我 们 过 去 用 来 计算 Kya 值 的 公式 。 

值得 注意 的 是 解 方 程 x** 一 a= 0 和 解 一 般 代数 方程 

ZoXx 十 ZIZx1 十 … 十 ax 一 0 

二 者 有 本 质 的 区 别 ， 对 于 方程 ** 一 a== 0 来 说 ， 初 始 近似 信 
*1 的 选择 是 无 关 重要 的 ， 不 管 选择 怎样 的 x*!， 经 过 若干 次 到 
近 计 算 以 后 ， 我 们 总 能 得 出 满足 精度 要 求 的 ka 的 值 。 而 对 
于 方程 (56 ) 来 说 ， 情 况 就 不 同 了 。 这 里 一 个 初始 值 可 能 得 
出 一 个 根 ， 另 一 个 初始 值 又 可 能 得 出 另 一 个 不 同 的 根 。 而 有 
些 初 始 值 根本 不 能 得 出 任何 确切 的 值 ， 也 就 是 由 公式 《62 ) 
计算 出 的 数列 x1,，x，，…，*m，… 不 趋向 于 任何 确定 的 极 
限 ， 即 是 发 散 的 。 


Xenk+1 一 %n 一 
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15. 微 商 的 几何 意义 


到 目前 为 止 ; 我 们 只 给 出 了 解 代数 方程 的 牛顿 方法 。 为 
了 把 这 种 方法 推广 应 用 到 解 任意 的 方程 上 ， 必 须 先 推广 微 容 
的 概念 ， 并 介绍 各 类 函数 的 微 商 。 为 此 ， 我 们 先 解释 微 商 的 
几何 意义 。 
考虑 多 项 式 
y=a0x* 十 a1xt-! 十 十 gx : 
的 图 像 。 任 取 两 点 KW 和 N 图 1 ) ，- 点 的 模 坐 标 是 z 点 


, 图 12 

的 模 坐 标 是 z+a。 则 点 MM 和 入 的 纵 坐 标 可 由 表达 式 

f(x%)==aox* 十 qix*-! 十 … 十 az 和 ?x 十 a) 二 as (x 十 a)* 十 a (4 

十 a)*-! 十 … 十 ax 分 别 计算 出 。 ee 一 割 线 并 计算 

它 的 斜率 Feo “从 图 上 可 以 看 出 t8y 一 全 ， 但 是 线段 MT 
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的 长 度 等 于 点 W 和 点 N 模 坐标 值 之 差 ， 所 MT 
mm Qo 
角 段 TN 等 于 这 两 点 众生 标 值 之 关 ， 所 以 TN 二 f(x 网 = fCz) 
由 此 得 出 
t sy 一 了 
Cs 
f(x+a)=f(x)+af’(x)+.…: 
其 中 ， 省 略 符号 表示 的 项 包含 a*，a’， ne 所 以 
tey= tt) + 


其 中 ， 省 略 符号 表示 的 项 包含 &，a?*……。 
因此 ， 割 线 MN 的 斜率 可 由 公式 
0 Ca) te “ (68) 


f(x+a)—{(a) 


& 


表示 。 po 
现在 ， 我 们 减 小 q 值 。 这 样 一 来 ， 割 线 将 绕 点 M 转 动 。 
' 当 a= 0 时 ， 就 是 极限 情况 。 此 财 割 线 成 为 曲线 4 二 f(x) 在 总 


图 13 


* 生 有 的 斜率 是 指 轴 肛 与 x 负 正方 向 天 名 角 的 正切 。 -例如 ， 加 果 二 条. 信 呈 
与 x 轴 的 亿 锋 度 为 60"， 则 它 的 余 率 等 于 tg60"， 即 v 3。 
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的 切线 。 图 13 中 给 出 了 a 一 1， 学 ， 二 时 的 情况 。、 


然而 当 a= 0 时 ， 公 式 《66 ) 中 点 省 略 的 项 都 为 零 。 所 
以 ， 多 项 式 g 一 Xz) 图 像 在 横 坐 标 为 点 的 切线 斜率 用 公式 表 
示 是 

ke=f (x) 
”因此 ， 多 珊 式 KKz) 的 微 等 于 多 项 式 的 图 多 在 模 人 村 为 

2X 点 的 切线 斜率 。 

例 . 求 多 项 式 1(x)==x* 一 4x? 十 5x 十 1 图 像 在 点 x 二 2 的 
切线 与 x 轴 的 夹 角 。 

因为 1 (x) 二 3x* 一 8x 十 5 

Ss (2)= 1， 所 以 tgy== 1， 即 夹 角 Y=45°。 
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“+6， 和 牛顿 法 的 几何 意义 


现在 ， 我 们 可 以 清楚 地 阅 明 牛顿 法 近似 解 代 数 方程 的 几 
何 意义 了 。 假 定 要 解 方程 人 (z)= 0 ( 人 xz) 为 多 项 式 ) ， 明 显 
地 ， 这 个 问题 就 是 求 函 数 ! 一 驴 z) 图 像 与 z 轴 的 交点 ， 即 是 求 
使 y 一 0 的 点 。 

假设 我 们 已 经 求 出 了 方程 根 的 近似 值 x,, 作 曲线 y 二 7(*》 
在 模 坐 标 为 zi 的 点 N 的 切线 ， 如 果 z: 选 择 恰当 ， 则 切线 与 * 
灿 的 交点 了 比 点 术 更 接近 函数 y 二 久 *) 与 * 轴 的 交点 ( 图 M)。 


图 14 
为 了 求 点 7 的 栅 誉 标 *,， 考 虑 三 角形 TMN， 三 角形 的 直 
角 边 MN 正好 就 是 函数 y==;(*) 在 点 *, 的 秆 ， 即 MN 二 f(x,)。 
另 一 直角 边 T 以 等 于 x 一 *，， 所 以 ， 切线 与 * 轴 夹 角 q ,的 正 
切 用 公式 表示 就 是 
tay = C68) 


和 %1 一 X2 
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人 公式 ( 68 ) 得 出 | 
(69) 


但 是 ，tgq91 是 曲线 yg 二 7(<) 的 图 像 在 横 毕 标 为 x; 的 点 的 切线 
。 射 康 ， 因 此 根据 微 商 的 几何 意义 ， tgp 一 f'(x1) 
所 以 ， 公 式 ( 69 ) 可 写 为 . 
i _ i(x1) 
Es en) 
于 是 我 们 就 求 出 了 方程 根 的 第 二 次 近似 值 z:,。 现在， 我 们 再 
作曲 线 ! 一 f(z) 在 模 坐标 为 z; 的 点 的 切线 ， 则 切线 与 ? 轴 交 点 
的 横 坐标 由 下 列 公式 得 出 - 
‘(x3) 
(xs) . 
、 一 般 地 ， 如 近似 值 *, 已 经 求 出 ， 为 了 求 下 一 次 近似 什 
as 4， 作 曲线 y 一 人 x) 在 机 前 标 为 x。 的 点 的 切线 ， 切 线 与 < 加 
交点 的 摸 坐 标 就 是 我 们 所 求 的 zsy\ 值 。 


xs 一 %z 一 


计算 x a+1 的 公式 是 : 
t= (70) 
或 这 
Anfi 一 za 一 区 ze (71) 
tgP» 


其 中 ,是 曲线 在 横 坐 标 为 zw 的 点 的 切线 与 轴 的 夹 角 。 这 个 


公式 与 牛顿 法 公式 ( 62 ) 完全 一 样 。 从 而 ， 我 们 已 经 知道 了 .， 


牛顿 法 的 几何 意义 、 它 实质 上 就 是 以 曲线 y= 人 *) 的 切线 代 、 
蔡 弧 线 。 由 于 这 种 原因 ， 和 牛顿 法 又 叫 切 线 法 。 | 


时 


图 15 中 说 明了 用 牛顿 法 得 出 的 点 列 x4， #3， "xo，…- 
是 怎样 表 近 曲线 yg 一 fx) 与 * 轴 的 交点 8。 


图 15 


179， 任 意 函 数 的 微 商 


前 面 对 和 牛顿 法 的 几何 解释 ， 能 够 推广 应 用 到 任何 具有 
#(*)== 0 形式 的 方程 ， 这 里 f(z) 可 以 是 多 项 式 以 外 的 函数 。 
为 求 这 种 方程 的 近似 解 ， 取 方程 的 根 的 某 个 近似 值 zx， 作 曲 

Me eh 

。 然 后 再 作曲 线 y 二 f(x*) 在 模 坐 标 为 x 的 点 的 切线 ,等 等 。 

pg yy 


pC?) 5 
YXn+1 一 和 mn te (71) 


其 中 tgqn 是 曲线 J 二 1(x) 在 模 坐 标 为 x 的 点 的 切线 斜率 。 
公式 ( 71 ) 仍然 不 能 用 来 进行 计算 ， 因 为 我 们 还 不 知道 
怎样 求 tgp»。 所 以 我 们 必须 先 学 会 计算 任意 函数 y=f(x) 图 
像 (不 仅仅 是 多 项 式 的 图 像 ) 的 切线 斜率 。 首 先 ， 我 们 先 计 
算 割 线 的 斜率 ， 过 函数 9 一 f(z) 上 一 点 M， 作 人 割 线 MN， 应 
用 前 面 多 项 式 的 论证 ， 我 们 同样 得 到 计算 割 线 斜率 的 公式 ， 
koi=tey= {e+e) f(s) (72) 
其 中 x 是 点 以 的 横 坐 标 值 ，x 十 a 是 点 N 的 横 坐标 值 。 如 果 减 
小 a 值 ， 割 线 将 绕 必 点 转动 ， 直 到 它 处 于 曲线 y==f(*) 在 这 点 
的 切线 位 置 ( 见 图 12 ) 。 所 以 我 们 可 以 把 切线 斜率 写 为 


Hm mate)-—Hs) 《73) | 


kis=tgPn= @ 


6 办 


s 


我 们 称 上 上 式 者 端的 极限 ; 为 函数 f(*) 的 微 商 ， 用 f'(x) 表 示 、 
即 


f(x)= Lm (+e) As) | | (74) 


现在 ， 我 们 可 以 把 (73 ) 式 写 为 
‘ks=tgpn=f (x) (75) 
| 所 以 任意 函数 ( 不 仅 是 多 项 式 ) 存 某 点 的 微 商 1'(x) 就 
是 曲线 ji 一 f(xz) 在 这 一 点 的 切线 斜率 "。 
因为 tg 二 f(xa)， 所 以 公式 (71 ) 可 写 为 
sat (76) 


a 
这 个 公式 与 公式 ( 62 ) 一 致 于 是 牛 顿 法 就 推广 应 用 到 了 所 
有 具有 f(x*)= 0 形式 的 方程 。 


bb 


所 如 果 不 能 作出 函数 了 一 f(x) 在 横生 标 为 x 的 点 的 切线 ( 比 如 函数 图 做 在 
访 一 点 有 间断 ) ， 则 革履 在 这 一 点 没有 福 商 。 - 


18. 微 商 的 计算 


、 前 一 节 里 ， 我 们 已 经 看 到 ， 为 了 求 与 曲线 9 二 帮 z) 相 狠 
的 切线 斜率 ， 必 须 计算 被 根 

3 f= tm zte) {(x) 

一 般 说 来 这 种 计算 非常 困难 。 但 是 ， 许 多 重要 实例 的 极 
限 已 经 计算 出 来 了 ， 换 名 话说， 就 是 最 常用 函数 的 微 商 我 们 
已 经 知道 了 。 下 而 列 出 了 一 些 常用 微 商 的 计算 公式 。 
(a)’=0 
(x2*)’ =kx*-! * 
(a7)’=arisa 
(sinax)’=acosax 
(cosax)’=—asinax 


@ QO@OO 


a 


ZJ) 一 一 一 
(tgox) COs*ax 
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a 
ctgaxz) "一 一 一 -一 一 
(ctgax) sin?*ax 


冯 


pe 
(los.z) zli@ 
(©) (arcsinax)! = Si 


a - 
(arctgax)’ = TT 
( 公式 @ 和 公式 @ 中 的 Ina 表 示 以 c= 二 2.71838… 为 鹿 的 对 数 ， 
称 为 自然 对 数 或 者 纳 皮 尔 对 数 ) 。 公 式 @ 中 的 不 仅仅 是 自 
然 数 ， 还 可 以 是 任何 实数 。 例 如 : 
1 .去 -1_ 


(V2) = ) 


CD a 


仅仅 用 公式 @ 一 @ 还 不 能 够 计算 出 所 有 函数 的 微 商 。 但 
是 ， 如 果 f(x) 是 由 一 些 已 经 知道 怎样 求 微 商 的 函数 经 过 算术 
运算 组 成 ， 则 我 们 还 是 可 以 很 容易 地 求 出 f(x) 的 微 商 。 为 
此 ， 我 们 使 用 下 面 的 法 则 ， 这 些 法 则 ( 和 公式 D~@ ) 在 高 
等 数学 课程 里 都 有 证 明 。 

@ 函 数 和 的 微 商 等 于 它们 微 商 的 和 。 

(F(x)+f x) = x)+f (2) 、 

加 常数 可 以 提 到 微 商 符号 外 。 

(af(*))’=af’(x) 

@ 两 个 函数 乘积 的 微 商 计算 公式 是 
E (fi(zjfs(z)) =7 x) (x) + x)’ (x) 

@ 分 数 函数 的 微 商 计算 公 式 是 

NONMEMGOMNOmNGOME 
f(x) Cf C2)) - 

在 第 13 节 中 给 出 的 多 项 式微 商 计算 公式 就 是 法 则 @ 和 法 

则 @ 以 及 常用 公式 @ 的 应 用 结果 。 


» _3x? 一 x 十 1 
例 1 求 函 数 人 (x) 一 2s 二 5 的 微 商 。 
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应 用 法 则 @@， 我 们 得 到 、 
fr(z)= (34 —#+1)’ (C20 a ee 
再 应 用 多 项 式微 商法 则 
(3x*:—zx++ 1 ) 一 6z 一 1 
(2xs 十 5) 一 6X: 
所 以 


f(x)= (2x3 + 8)(6i— 1) (g++ : 中 
(22 十 5)2 . 
一 二 6%“ -Fat 一 6x* 十 30x 一 5 
(22 45). 


例 2 求 本 数 Kx)= 直 (aresinsx 一 汉 ) 的 微 商 。 
应 用 公式 @ 和 公式 @， 以 及 法 则 中 和 @， 我 们 得 到 
(x)= 直 


[re -二 (到 也 ) 


3 
T1071 5 

例 3 求 函数 f(x)=107sin2x 的 微 商 。 

应 用 法 则 @ 和 公式 @、 多 得 ， 
f(x)=(107)’sin2x+10(sin2x*)’ 
=10xsin2x*In10 十 10z2cos2x 

一 10xz(sin2x*ln10 十 2cos2x)》 

上 面 介绍 的 求 微 商 的 法 则 ， 已 经 能 够 计算 出 许多 函数 的 
微 商 。 还 有 另外 一 种 非常 重要 的 法 则 一 一 计算 复合 函数 微 商 
的 法 则 ， 用 公式 描述 如 下 : 

?1 


如 果 函 数 4 一 所 2z) 可 以 写 为 3 一 P(z7》 z=q(t) 的 形式 。 
其 1(*) 的 微 商 计算 公式 是 
f(x2)=F'(:)"m’(*) (77 ) 
其 中 z 一 9(z7。 
例 准 泛 数 g 二 sin(x ) 购 牙 商 。 

”这 个 现 玄 可 写 为 二 sinz，: 二 Xs。 函数 F(;) 二 sin: 的 微 
商 是 F/(:) 一 cos:， 函 数 p(z) 一 2 的 徽 商 是 9'(2) 一 3z*。 由 
公式 (77 ) 得 : 

{sin(x*))’=F’(s:).9’ (x)=cosz*3x: 
将 z==** 代 入 ， 我 们 得 到 
(sin(%x’)) ‘=3x?cos(**) 
关于 微 商 的 概念 ， 读 者 可 在 其 它 书 中 找到 更 详细 的 说 
组 。 例 如 证 Ya.B , 泽 尔 多 维 奇 的 《 高 等 数学 入 门 》 一 书 里 就 
和 详细 说 明 。 


19. 求 初始 近似 值 


现在 ， 我 们 讨论 初始 近似 值 的 选择 问题 。 求 解 方 程 {(x} 
呈 0 时， 第 一 次 初始 近似 值 可 通过 图 示 法 得 出 ， 但 对 此 要 作 
函数 y 二 光 x) 的 图 像 ,并 求 出 它 与 * 轴 的 交点 《 在 交点 上 4 二 0， 
所 以 Kxz)=0 )。 站 


如 在 计算 机 - 饰 方 程 ) 。 因 此 ， 必 须 采 用 另外 的 求 初始 近似 
值 的 方法 。 鸭 方 法 通过 计算 函 激 在 某 些 点 -上 的 西数 值 


( 例如 让 和 变量 取 定义 范围 内 的 整数 值 ) 来 确定 初始 值 的 先 
择 。 如 采 函 数 y 一 并 z) 连 续 ( 即 函 数 图 像 没 有 间断 点 )， 则 ” 
当 琐 数 在 自 变 量 分 别 取 a 和 8 信 有 相 异 的 符号 时 ( 图 16a ) ， 
在 a 和 之 间 有 方程 (x)= 0 的 根 。 此 时 ，a 和 5 都 可 作为 方程 
f(x)= 0 根 的 初始 近似 值 。 如 果 函 歼 图 像 有 间断 点 ， 也 可 能 
出 现 函 煞 从 负 值 不 经 过 零 一 下 跳 到 正 值 的 情况 ( 图 165 ) 。 


上 2 1 
(a) (DD (0 


要 注意 的 是 这 种 方法 可 能 遗漏 掉 方程 的 某 些 根 。 例 如 图 
3 


16c 中 出 现 的 情况 ， 函 数 g=1(*) 在 a 和 4 两 点 有 相同 的 符号 
但 是 在 a 和 b 之 间 却 有 使 方程 为 零 的 点 。 

于 是 ,我 们 已 经 求 得 两 点 :4 和 ?。 应 该 选 娜 一 点 作为 牛 
顿 法 的 初始 近似 值 z, 呢 ? 考虑 图 17， 图 17a 和 图 175 表 明 ， 恕 
果 函 数 曲线 是 向 上 目的 ， 就 选择 使 函数 为 正 值 的 点 ，4 点 或 8 
点 作为 初始 近似 值 。 初 始 值 选 择 不 当 就 可 能 导 儿 计算 出 的 z。 
值 超出 区 间 (a, .5 之 外 。 同 样 地 ， 如 果 函 数 笛 线 是 向 下 站 
的 (图 17c ,174) ,就 选择 使 函数 为 负 值 的 点 作为 初始 近似 值 。 


(b) 


7d) 


图 17 


如 果 知 道 函 数 3 一 人 (xz) 的 图 像 ， 使 用 这 条 规 则 是 很 方便 
的 。 但 如 果 没 有 作出 函数 的 图 像 ， 则 需要 计算 函数 {(x) 的 二 
阶 微 商 来 判断 函数 的 凹凸 性 。 画 数 xx ) 的 二 阶 微 商 就 是 它 一 
阶 微 商 的 微 商 。 
例如 ， 给 定 函 数 江 zx) 一 2 一 422 十 3x 一 1， 它 的 一 阶 微 
商 是 f'(*) 二 3x: 一 8 十 3， 而 二 阶 微 商 是 f(x*)=6x 一 8。 
74 一 


了 


在 区 间 (a，5) 上 ， 如 果 函 数 的 二 阶 微 商 是 正 的 ， 则 还 
数 在 比 区 间 上 是 向 上 媚 的 。 反 之 ， 如 果 二 阶 微 商 是 负 的 ， 财 
函数 在 此 区 间 上 是 向 下 目的。 这 些 在 高 等 数学 课程 中 都 是 有 
证 明 的 。 运 用 这 些 知识 ， 我 们 得 出 下 面 使 用 牛顿 法 的 规则 

设 函 数 i(x) 在 区 间 端 点 a 和 48 有 相 异 的 符号 ， 如 果 信 x) 在 
区 间 [oe，8】 上 的 二 阶 微 商 是 正 的 ， 则 应 选取 使 函数 从 *) 为 
正 值 的 点 ，a 点 或 5 点 作为 初始 近似 值 z:。 反 之 .如果 fx) 在 
区 间 〔o，"] 上 的 二 阶 微 商 是 负 的 ， 则 应 选取 使 函数 儿 x) 为 
负 值 的 点 作为 初始 近似 值 zke 


人 


20. 人 法 


“在 求解 方程 中 ， 和 牛顿 法 常 结合 着 纺 截 法 一 起 使 用 。 如 
王 数 3 一 f(z) 的 图 像 是 向 上 外 的 ， 则 求 o 和 zi 分 别 用 公式 - 


Qi 一 0 一 fo) = ta) (78) 
ID pa 
x1=b PC) C29) 


反之 ， 如 果 函 数 #=1(*) 的 图 像 是 向 下 四 的 ， 则 求 a 用 公式 
《78 ) ， 求 x 用 公式 
f(a) 
f'(a) 

如 图 18a 和 i185 中 所 见 到 的 那样 ， 方 程 a f(x) = 0 的 根 & 通 
常 位 于 点 ci 和 xz: 之 间 。 然 后 再 应 用 牛顿 法 和 弦 截 法 ， 我 们 又 
得 出 一 组 新 的 点 6: 和 za， 如 此 等 等 。 


Xi 二 4 一 


(80) 
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下 这 种 方法 得 到 两 组 点 列 : a1，a;，…，cr，…: 和 #1y 
zx:，…，xa，…。 这 两 组 点 列 从 不 同 的 方向 逼近 方程 的 根 8 
这 种 方法 的 好 处 就 是 在 于 它 分 别 从 根 的 两 边 -- 起 逼近 。 ”. 

例 用 组 合法 解 方程 x 一 sinx 一 0.5= 0， 再 求 精度 为 
0.001。 | 

作 连 续 函 数 {(x) 二 x 一 sinx 一 0. 5 的 后 数 值 表 : 


由 表 可 知 ， 方 程 的 a a rE 
根 在 1 和 2 之 间 。 应 用 |- 0 2 


第 18 节 中 的 公式 2 和 公 f(x) Ca |-o. | —0, 5 上 
式 4 得 ff (zx) 一 1 一 cosz， 
产 以 在 这 里 牛顿 公式 为 
2 一 人 一 Sinza 一 0.5 ( 81) 
1—~cosxn 
为 了 确定 应 选取 1 或 者 2 作为 初始 近似值 xs， 我 们 求 函数 
fx) 的 二 阶 微 商 。 应 用 第 18 节 中 公式 @@， 它 的 二 阶 微 商 是 
人 (x) 二 sin*， 且 在 区 间 (1，2 ] 上 ，sinx 是 正 的 "。 明 显 
地 ， 应 用 前 面 叙 述 的 规则 ， 应 该 选取 2 作为 初始 近似 信 x。， 
因为 Kx) 在 点 2 的 值 是 正 的 。 
应 用 公式 (81 ) 有 


_，_2 一 sin2 一 0.5- , .2—0.909—0.5. 
mn re 
=1.583 
昂 一 方面 ， 应 用 公式 {《 78 ) ， 我 们 又 得 到 


二 工 一 (一 I ed 
ai 1 —(~0.341)0 69100-311)™1.366 


* 因为 sinx 在 区 间 【0, x = 〔0，3.141…] 上 是 正 的 。 所 以 sinx 在 
它 的 于 区 间 1 ，2 上 也 是 正 的 。- 
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然后 对 区 间 Cas, | ii ( 81 ) 和 公式 (78) 
我 们 得 到 
xa m1.588-.1.583—1.000—0.5 


T0012 一 1.501 


1.583 一 1.366 _ 
0. 083+0. 113 


2a 一 1.366 十 0.113 
继续 下 去 ， 我 们 得 出 
xX3=1.498 
Cs 一 1.498 、 

因此 ， 方 程 精确 到 0.001 的 根 是 1.498。 


一 1.493 
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适 代 法 的 收 伍 检 验 


现在 我 们 应 用 微 商 的 概念 ， 引 出 一 种 新 的 选 代 法 收敛 检 
验 。 为 此 ， 需 要 拉 格 朗 日 ( Lagrange ) 公式 ( Lagrange 是 
18 世 纪 法 国 数学 家 ) 。 

考虑 区 间 【a，b] 上 曲线 g 一 Xz)， 设 M 和 六 分 天 表示 
曲线 在 此 区 间 上 的 起 点 和 终点 ， 作 弦 MN ( 图 19 ) ， 则 弦 的 
斜率 是 

上 弦 二 tg 中 二 这 


但 是 ，MP=b 一 a， PN= f(5) 一 f(a)。 所 以 
f(b)—f(u) 


飞 到 二 b—a 


图 19 
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用 T 表 示 从 弦 MN 到 弧 MN 上 距离 最 远 的 点 ,过 这 一 点 作 : 
弦 MN 的 平行 线 ， 它 必 将 与 曲线 相 切 。 因 为 如 果 它 与 曲线 栓 
交 的 话 ， 就 说 明 在 曲线 上 还 有 另外 的 点 茅 芝 WMA 的 距 离 比 点 
T 到 弦 WN 的 距离 远 。 换 句 活 说 ， 就 是 曲线 在 了 点 的 切线 平行 
于 弦 MN， 因 此 它 与 弦 有 相同 的 斜率 。 然 面 ， 切 线 的 入 率 等 
于 #(c)。 这 里 < 是 点 7 的 席 坐 标 ， 所 以 ， 我 们 得 出 公式 


fc) 一 让 的 二 区 人 (82) - 


这 个 公式 就 叫 拉 格 朗 日 公式 。 要 注意 的 是 拉 格 朗 日 公式 
中 ， 点 "总 是 介 于 点 ac 和 4 之 问 的 。 拉 格 朗 日 公式 也 可 写成 下 
面 的 形式 

f(6)—f(a)=f’(c)(b—a) (83) 

现在 ， 我 们 再 回 到 用 迭代 法 解 方 程 +=qp(%) 上 来 。 假 定 
映射 y= 二 q(x*) 把 区 间 Ca，&) 映射 到 它 自己 的 子 区 间 ， 则 在 
这 区 间 上 有 不 等 式 | 9 (xz) | 之 4， 其 中 4<< 1 。 因 而 在 区 间 
[a,b】 上任 取 两 点 + 4 和 xs， 点 9(x1) 和 q(x) 也 属 于 区间 
[a，b】， 应 用 拉 格 朗 日 公式 下 

q(x2)—9(x1)=9 (cs —x,) 
其 中 c 介 于 x1 和 x 之 间 也 属于 区 间 (g,， 85), : 且 | 9 (ce) | 过 
&< 1 成 立 。 所 以 

| 9(*2)— 9(*1) | 4 1x2—x1 | (84) 

不 等 式 ( 84 ) 说 明 ， 9(z) 是 一 种 压缩 映 射 。 但 我 们 知 
道 ， 如 果 x 一 >qp(*) 是 把 区 间 { a，64 映射 到 它 自 已 子 区 条 
的 压缩 映射 ， 则 对 区 间 上 任意 点 *。， 点 列 zo，X1，Xa， 
Xa 《其 中 xn+1 二 (xn) ) 收敛 于 方 织 xz 一 (xz) 的 根 。 于 是 
我 们 就 证 明了 下 面 的 定理 。 

定理 设 函 数 y 二 q(x) 把 区 间 (ae，4] 了 瑞 射 到 它 自 己 舱 
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子 区 间 。 且 在 区 闻 上 育 不 等 式 | 9(z)| <s(4<1) 成 立 。 
则 对 区 间 上 任意 点 ze， 由 zari 一 (zs) 得 出 的 点 列 x*。，*1 
… Xn， … 有 收敛 于 方程 ?一 9(z) 的 根 。 
一 往 咯 地 讲 ， 刚 才 证 明 的 定理 说 明 ， 逐 次 避 近 法 使 我 们 能 
够 求 出 方程 :二 p(*) 的 那些 满足 不 等 式 | 9'(&) | <1 的 根 8 。 
法 代 法 的 几何 描述 ( 见 第 8 节 ) 就 是 这 些 点 连 成 一 条 折线 ， 
当 19 (8) | > 1 时 ， 就 破坏 了 这 条 折线 。- 
如 果 在 整个 数 轴 . 上 满足 不 等 式 ,| p(x) 1 <d< 1， 则 
选 代 法 的 收敛 不 依赖 初始 近似 伪 ze 的 选择 人 《 见 第 10 节 ) 。 
例 1 兴 代 法 能 否 用 于 解 下 方程 ? 


-29SY 直 sin 


:4 
这 里 9 一 casx 二 Sing， 所 以 
“(x)= sinx+coss 
| p(x)= 4 
得 -| sinz | 和 1，j| cosx | 寺 1， 因 此 3 _ 
lq9’(Cx)| =|=ns+ cosx| lsin 旭 十 lcos 双 
了 4 ' 4 


1 
.<3 
产 凡 能 用 潜 代 法 解 方程 2 一 9 时 地 于。 
例 2 ”能 否 用 迭代 法 解 方程 ( 85 ) 
xz 一 在 一 2 (85) 


所 求 很 在 区 闻 ( 1 ， 2 上， 因为 连续 函数 4 一 :一 4 十 3 
在 这 个 区 间 上 改变 符号 ， 即 1 一 4 十 324< 0，2 一 奸 2:>>0 儿 
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这 里 p/(z)= 一 2*1n2。 我 们 再 计算 表达 式 一 2*In2 在 区间 
[1，2) 上 的 值 。 如 果 1<#<2， 则 2 所 27 人 4。 所 以 
21n2<2*1n2::. 41n2。 . 
从 自然 对 数 表 (底数 是 ce 一 2.78 ) 中 查 出 tn2 一 0.69…， 最 
然 在 区 间 [a，b】 上 ， 不 等 式 1.38 委 2zin2<<2.76 成 立 。 所 
-以 不 能 使 用 选 代 法 。 
为 了 能 够 使 用 选 代 法 ， 我 们 把 方程 (85 . 写 为 
2 一 4 一 2 
的 形式 ， 然 后 两 边 取 以 2 为 廊 的 对 数 ， 则 我 们 得 到 
Y 一 Jog:(4 一 2z) 
这 里 
Ue oR 
9 2) 一 (4 一 z)in2 
县 在 区 间 (1，2 ] 上， 不 等 式 


19 (xz) | < -=< 1 


下 1. 和 


威 立 〈 读 者 可 很 容易 地 导出 这 个 不 等 式 ) 。 
“ 所以， 方程 写成 这 种 形式 时 ， 使 用 夺 氏 法 收敛 。 
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22， 适 代 法 的 收 伍 速 度 。 


现在 ， 我 们 利用 王 数 9(*) 的 微 商 求 先 代 法 解 力 入 
x 二 q(x) 的 收敛 速度 。 即 我 们 要 求 方程 根 E 的 近似 值 x ,, x，， 
…，zn， 与 8 的 误 益 ww 一 一 Zn 减少 的 速度 。 
必须 注意 到 方程 8 一 9(E) 和 Xn+ :站 (zo) 部 是 准确 的 。 
由 此 有 下 面 的 等 式 
-aanrm 一 5 一 znrlm9(E) 一 (zw) 
但 由 拉 格 朗 日 公式 有 
9(E) 一 9(zn) 一 9 (cn)(E 一 zw) 一 9 (ca)an 
其 中 c* 介 于 xn 和 5 之 间 。 所 以 
Qnr1 一 申 - ‘(cn)an 《86 ) 
从 等 式 ( 86 ) 可 引出 下 面 的 结论 : 
设 E 是 方程 x 二 p(x) 在 区 间 (a，b) 上 的 根 ， 如 果 在 这 个 
及 间 上 ， 不 等 式 | 9'(x) | 之 4 之 1 成立， 县 初始 近似 值 * 
也 属于 区 间 (a，8) ， 财 对 任何 x 有 关系 式 | 


| acer i<a* la 1 3 (87) 
成 立 。 . 
实际 上 ， 由 等 式 (86 ) 知 
las l= 19’Ce)1: 1ail 


但 c: 属 于 区 间 [a，# )《 图 20 ) ， 所 以 


站 这 一 节 对 切 学 者 可 以 省 路 。 


leeD 1<a 


由 此 得 出 
Ja: ij<ea ic | 
a 
图 20 


` 用 同样 的 方法 ， 我 们 可 得 到 
{asl=Ip ‘(clasl<slasl<a lal 
一 般 地 
| anf | <a" las | 
这 就 证 明了 我 们 的 论述 。 
因为 0<&< 1， 随 着 4 增 大 ， 数 列 4，4*，…，4*… 六 
向 于 零 ， 误 差 svri 也 这 庙 于 替 。 换 名 话说， 上 面 假 定 的 二 
似 信 x,，x,，…，%。… 通 近 时 其 误差 | 一 za+1 | 减少 
比 | ak | "人 块 。 
同样 我 们 可 以 证 明 ， 如果 在 区 间 ? o, 8 ) 上 ， 不 等 式 
jp'(x) | > 1 成立 。 则 和 迭代 法 发 散 。 
- ”如 果子 数 的 微 商 8“(z) 在 E 处 为 惟 ， 则 迭代 法 收敛 的 速 
度 最 快 。 这 样 ， 随 车 我 们 得 到 的 值 接 近 上 ，qp'(x) 就 趋向 于 
零 。 因 为 
(am 9 cov ion | 
所 以 ， 和 迭代 法 收敛 的 速度 随 着 通 近 8 而 增 快 。 
我 们 用 鞍 代 法 求 平方 根 时 ， 已 经 过 到 过 类 似 的 情况 。 记 


得当 时 我 们 是 用 方程 zs= 工 :2 代替 方程 za=a。 然而 函数 
9(z) 一 志士 2 的 微 商 是 


(Xa) 2x— (tta)(2x)! 
9/(z) 一 < 
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_2x。2z 一 (23 十 ao) 2_ zz 一 2 


4X2 2x7 
《 见 第 18 节 规划 @ 和 第 13 节 公式 (55) ) ， 所 以 
a):a_ 0 
vy'(vVa)= 2 4~0 


因此 ， 函数 p(x) 的 微 商 在 x 二 VG 点 时 为 零 ， 并 且 随 著 逼近 
点 zx 一 va 时 ， 加 快 了 方法 的 收敛 。 

- 随 着 禹 近 方 程 的 根 而 加 快 收敛 速度 ， 这 一 点 也 是 牛顿 法 
的 一 种 特性 ( .上面 求 平方 根 的 方法 就 是 一 种 特殊 实例 ) 。 实 
际 上 ， 我 们 知道 : 牛顿 方法 就 是 用 方程 


代替 方程 Kz)= 0， 稀 局 各 用 寻 次 过 放 解 演 个 力 生 。， 于 是 
我 们 有 
f(x) 
| 
而 
9(s)=1— [ 作 的 ] 


二 ee 
了 4 x) 2 


(f(z) ’:—f(x)f"(*) 
Fx 
f(x)f"(x) 
CF C2))? 


镶 为 方程 在 点 上 满足 人 E) 一 0， 故 9/ (8) 一 0。 所 以 ， 如 前 


面 所 述 的 那样 ， 获 着 逼近 点 8 而 加 快 了 方法 的 收敛 速度 。 
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23， 逐次 从 通 近 法 解 线性 方程 组 


到 现在 为 止 ， 我 们 已 经 学 会 了 解 一 个 未 知 量 的 方程 。 下 
面 我 们 将 讨论 一 次 方程 组 的 解法 。 
假设 有 mx 个 未 知 量 的 mm 个 一 次 方程 *: 


QliXit 十 G12%2 十 十 dimXm 二 1 
[> 41X1 十 az2X2 十 … 十 azmXm 一 bx . (88) 


QmiX1T ans x 二 十 QammXm== bm 


这 样 的 方程 组 有 非常 广泛 的 应 用 ,例如 ,测量 学 工作 者 综合 地 
球 表面 大 面积 改造 的 测量 结果 ， 他 们 有 时 必须 解 几 百 个 方程 
的 方程 组 ， 工 程 师 们 设计 刚性 结构 以 及 其 他 专家 们 在 许多 特 
殊 的 领域 中 ， 也 必须 解 这 样 的 方程 组 。 
解 这 类 方程 组 常用 的 方法 ( 如 未 知 量 消去 法 ) 一 般 都 非 
常 及 烦 ， 使 用 逐次 逼近 法 求解 就 比较 方便 。 首先 我 们 看 一 个 
用 逐次 焉 近 法 解 方程 组 的 例子 。 
假定 有 方程 组 
10x1 一 2x, 十 x*s 二 9 
| Xi 十 5z， 一 zs 一 8 
4x1 十 2x: 十 8xs 一 32 


上 * 在 这 类 方程 里 ， 我 们 用 字母 xi，x2，.…，xw 交 示 米 知 量 ， 用 字 母 a 1 3 


网 示 系数 ， 其 中 系数 的 第 一 个 下 标 表示 方程 的 个 数 ， 第 二 个 下 标 表 示 未 知 量 的 
个 数 。 例 如 ak7 均 示 第 四 个 方程 第 七 个 变量 的 系数 。 
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要 求 出 未 知 量 z:、 za 和 2%s 的 值 ， 精度 为 0.01。 

我 们 从 方程 组 的 第 一 个 方程 中 解 出 z:， 从 第 二 个 方程 中 
解 出 x,， 从 第 三 个 方程 中 解 出 x*s， 则 方程 组 为 ， 

| [0.9+ 0.2x: 一 0.1xs 

x 二 1.6 一 0.2x1 十 0.2xs (89) 
一 4 一 0.52: 一 0.25z: 

任 取 一 组 全 作为 x,、*s 和 xs 的 初始 近似 值 ,比如 取 x,'"*== 0， 
x2' 二 0，xs'"} 二 0， 把 这 组 值 代入 方程 组 (89 ) 的 右 端 
并 计算 可 得 到 xz:、x: 和 xs 的 下 一 次 近似 亿 ， 


2 一 0.9 

(a : 
xa‘!1?== 4 

把 得 出 的 这 组 近似 值 再 代入 方程 组 (8 89 ) 右 端 计算 又 得 到 下 


一 次 近似 值 : 
OB te Btn 4 一 0.82 
X22 =1.6—0.2xX0.9+0.2X 4 一 2.22 
xa‘ =4—0.5x0.9—0.25x1.6=3.18 
一 般 地 ， 如 果 近 似 值 x1‘”?，x，' "和 x “" 已 经 求 出 ， 则 要 
求 下 一 次 近似 值 就 使 用 公式 - 
ne .9 十 0. 2x:5n) 一 0. 1zs°”2 a 
x "t=1.6—0.2x1.™ +0.2xs." (90) 


xa‘ "+13~ A —0.5x1 "7—0.25x2."? 


此 例 中 各 次 计算 结果 如 表 2 。 
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oe | 2 :， 3 | 4 | 5 | 6 
| | 0,9 0.82 | 1.03 | 1.01 Ti 00 
an | 1.6 | 2.22 or | 2.00 | 1.99 | 2.09 

A 5 - | 二 | 
xa(1) | 4.0 | 3.15 | 也 03 | 2.07 、| 3,00 | 3.00 

| 


我 们 看 到 ， 在 精度 要 求 内 ， 下 列 等 式 成 立 ， 
(1 x 


(i 2 (91) 
Xa 558 一 和 6869 


在 方程 组 ( 90 ) 中 下 w= 5 并 考虑 等 式 (91 )， 在 精度 要 求 
内 我 们 得 到  - 
[和 1 50.9+0.2x2°5° —0,1x:5’ 
和 55 1.6 一 一 0.22) 人 0 32 
zs' "v4—0 5x1553 一 0. 25z:55) 


( 事实 上 ， 等 式 是 精确 的 ， 不 过 这 并 不 要 紧 。 ) 由 此 得 出 
X11.00，x*a，， 王 2.00，xs'， 二 3.00( 在 精度 要 求 内 ) 
煽 方程 组 的 握 Fe 

一 髓 的 方 各 组 也 可 用 同样 的 方法 求解 "。 候 定 有 方 程 组 
( 88 ) ， 从 第 一 个 方程 中 解 出 *,， 从 第 二 个 方程 中 解 出 xs。' 
如 此 等 等 。 则 ( 88 ) 式 为 


> 这 一 节 的 以 下 部 分 ， 对 于 初学 者 可 以 省 略 。 
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f a G1 
A 2 
1 Gii Qs QI 
! 和 Gam 
| za 一 二 人 
/ G22 G22 - G22 (92) 
| 
了 A 未 
x RY 
tmm Gmm Umm 


设 zz，…，xorcl3 为 未 知 量 4a，xa，…，zm 的 一 
组 初 纵 近似 人 把 这 组 值 代入 ( 92 ) 式 右 端 计算 ， 我 们 得 到 . 


所 求 未 知 量 的 第 二 次 近似 值 

人 各 一 aa 
Qari g QI11 LE 
b 0 

jz 一 7 222171513 一 一 一 
Q22 022 - 03 

| 用 

a a - 

(nsa9 一 -一 miy 1 — so 

Qnm Gmm Qmm 


同样 地 ， 如 果 第 x 次 近似 信 已 经 求 出 ， 则 下 一 次 近似 值 的 计 
。 算 公式 是 


ptt bt 1 — lg 
Qi Q11 QIL 
恩 asi _a 

Xr tl 一 21X on3 一 gn 2m ony 


Qs2 022 Q22 
必 a Qmsm-— 
| mi "(93) 
Gmm Qnm amm 人 


现在 ， 我 们 考虑 一 个 与 逐次 逼近 法 解 方程 组 光 关 的 问 
题 ， 即 是 用 解 方程 组 的 方法 求 某 种 逼近 过 程 的 极限 状态 。 


设 有 3 只 水 桶 ， 第 一 只 桶 里 装 有 12 升 水 ， 第 二 只 和 第 三 
只 是 空 柄 。 先 把 第 一 只 桶 里 的 水 倒 一 半 到 第 二 只 桶 里 ， 然 后 
把 第 二 只 桶 里 的 水 倒 一 半 到 第 三 只 桶 里 ， 最 后 再 把 第 三 只 桶 
重 的 水 倒 一 半 到 第 一 只 桶 里 。 这 样 重复 循环 倾倒 20 次 后 ， 求 
每 只 桶 里 各 有 多 少 公升 水 ( 精度 为 0.0001 升 ) ? 

明显 地 ， 这 个 问题 就 是 讨论 逐次 通 近 水 分 布 的 极限 状 
态 。 这 种 极限 状态 表示 循环 倾倒 水 不 再 改变 水 的 分 布 情况 。- 
如 果 开始 循环 倾倒 时 ， 第 一 只 桶 里 有 xz 升水 ， 第 二 只 桶 里 有 3 
”升水 ， 则 第 三 只 桶 里 就 有 ( 2 一 x 一 9 ) 升水 ( 水 的 总 数 在 候 
例 过 程 中 不 改变 ) 。 则 循环 一 轮 倾倒 的 情况 如 下 ， 

第 一 只 桶 ， 第 二 只 桶 。 第 三 贝 桶 


初始 状态 x y 12 一 ZX 一 四 
第 一 次 倾倒 后 由 冯 +3 12 一 zx 一 ?9 
rg . 2 XY ot Pi 
第 二 次 倾倒 后 3 4 + 12 本 2 一 总 
一 - 人 和 en 
第 三 次 倾倒 后 6 十 各 一 各 。 到 十 三 6 一 吕 z 一 上 
如 果 这 样 倾倒 后 ， 每 只 桶 里 的 水 数量 不 变 ， 即 湾 足 下 面 
方程 C : 
= A 
| zx 二 6 十 5 一 二 
要 用 外 


解 此 方程 ， 我 们 得 到 x 二 6，y 二 3。 因 此 ， 水 的 极限 分 布 状 
态 是 第 一 只 桶 里 有 6 升水 ， 第 二 只 桶 里 有 3 升水 ， 第 三 只 桶 
里 也 有 3 升水 。 A 国 、 
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让 我 们 再 看 水 的 分 布吉 近 极限 状态 的 速度 是 多 少 。 假 设 
第 一 只 桶 里 有 a 升水， 第 二 只 有 5 开水， 经 过 一 轮 循环 侧 全 


. 后 ， 第 一 只 桶 里 的 水 是 
QQ (94) 


b (95) 
- 3， 则 


令 a=a 一 6，8=b 一 3，a:=a: 一 6， Bi 一 5 一 
由 等 式 (94) 和 (95 ) 得 : 


所 以 ， 如 果 |a | <e，181<e， 风 


13 、 、 
| a: |< 芒 =s0.2e 


18: | < 二 ss 0.34e 


这 说 明 经 过 两 轮 循环 倾倒 后 ,误差 4 和 8B 至少 减 小 为 1/3 僧 ， 
因此 经 过 20 轮 循环 倾倒 后 ， 误 差 至 少 减 小 到 1717000(35 
心 70000 ) 。 所 以 ， 衔 环 倾倒 20 轮 后 ， 在 精度 0.0001 内 水 的 
分 布 状态 是 ， 第 一 只 桶 里 有 6 升水 ， 第 二 只 和 第 三 只 桶 里 各 
有 3 升水 。 


24， 逐 次 逼近 法 解 非 线 性 方程 组 


隶 次 逼近 法 ( 迭代 法 ) 还 可 以 用 来 解 某 些 非 线性 方程 
组 。 例 如 ， 考 虑 方程 组 


x+y? (96) 


选取 x。 三 0，yo = 二 0 作为 初始 近似 值 。 把 这 组 值 代入 方 
程 组 等 式 右 端 计算 得 到 下 一 组 近似 值 *, 二 2，y1 = 二 1。 再 将 
这 组 全 代入 ( 96 ) 式 右 端 计算 ， 我 们 又 得 到 
[= - 士 L=2.25 


l= 1 +2 直 + =1.15 


继续 下 去 ， 我 们 可 得 出 
2 二 2.25- 二 .15 一 2.31 


本 2.25 十 1.15: 
3 一 1 0 


3 
LS 2 二 231 十 1.18 动 二 理 =2.33 


2.31 十 1.18: 
让 20 


一 1.18 


n=1 一 1.18 
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[2 337 Ee 18_, 33 


2. -33 18° 


[3 =1.18 


~“ 可 以 看 出 ， xX4=%6=2.33, ee 
所 以 方程 组 的 解 是 
pn 
y=1.18 
一 般 地 ， 如 果 友 方程 组 
人 3 ) 
b=Y(x, £) 
其 中 q(x，5 ) 和 u(*，2 ) 是 基 种 函数 。 我 们 选取 涪 始 近似 
值 *。 和 yo。 代入 《97 ) 式 右 端 并 按 公 式 
Xn+1 一 一 (Zi + sn ) 


(97) 


(98) 


fnt1=(xns Sn) 
F 茜 个 n 值 ,等 式 %nt i 从 Xny arm 
ln 4 一 fn 


2 pe 4) 是 定义 起 XX、 平面 


上 一 个 有 界 闭 # 数 ， 也 就 是 说 ， 假 设 区 域 D 全 部 在 


某 矩 形 的 内 部 并 二 办 点 。 比 如 图 、 乡 边 形 、 椭 圆 等 


此 外 ， 我 们 还 假设 函数 o(x，3# )， 
外 (7， 3 和 中 (2， 4) 在 区 域 D 上 定 
义 了 一 种 映射 ， 它 把 区 域 5 当 射 到 同一 平面 的 其 个 区 域 上 。 
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为 了 求 区 域 D 上 点 4o(*e，3o) 的 了 映射 点 将 x， J 值 代 入 
q(x，9) 和 中 (x*，#) 计 算 ， 其 结果 就 是 Me 点 的 象 的 坐标 。 例 
如 ， 假 设 
pd y) 一 2 十 可 
Vx, #) 一 2x3 
则 点 Mo( 1，3 ) 被 映射 到 点 No(10，6 )。 
以 后 ， 我 们 用 字母 表示 函数 q(x*，3) 和 只 x*，) 给 出 
的 映射 ，@(M) 表 示 点 M 在 映射 下 的 象 ，@(D) 表示 区 域 D 
作为 一 个 整体 在 映射 下 的 象 。 
假定 映射 把 区 域 D 映 射 成 它 的 子 区 域 D:= @(D)， 则 
应 用 同样 的 映射 ， 可 以 把 区 域 DD, 映 射 到 它 的 子 区 域 D,= 
中 (Di) 上 。 当 然 忆 :也 含 于 区 域 D 里 面 。 继 续 这 种 过 程 ， 我 
们 得 到 一 系列 一 个 含 于 一 个 的 区 域 D，D1，Ds，…，Dns* 
( 辐 21 ) 。 
* 如 果 对 于 区 域 卫 内 任 意 两 点 2 和 M,, 存 在 一 数 a 
(0<<i<< 1 ) ， 使 得 不 等 式 
;{ PM), DM,)) <i? (Ms 加) 
成 立 。 我 们 就 称 映射 中 是 一 种 
压缩 。 这 里 1(M，NN ) 表 示 点 
AM、N 之 间 的 距离 。 
就 单 变量 来 说 ， 下 面 的 命 
题 是 能 够 被 证 明 的 : 
假设 映射 中 把 区 域 石 映 射 , 
到 它 的 子 区 域 ， 且 @ 是 一 种 压 图 21 
缩 ， 则 在 区 域 了 上 存在 唯一 的 点 N，= AN @(N )。 并 且 这 点 属 
_ 于 所 有 区 证 D,， 点 的 坐标 5、? 潢 足 方程 组 ( 97 ) ， 即 
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7) 
1 一 由 85， 7) 
若是 单 变量 ，E 和 mn 就 是 由 迭代 法 计算 出 的 近似 值 。 如 果 


Mo(ze，3o) 是 区 域 D 上 任 一 点 ， 且 Ma+1= @( 即 x*a+1 二 9 
(zr, 3a)，3gn+1 一 由 xn，ga) ) ， 则 点 列 Mo，Mi Mu 
… 收 敛 于 一 固定 映射 点 N。 
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25， 修 正 距 离 


迭代 过 程 的 收敛 性 ， 要 求 是 一 个 压缩 映射 仅 是 一 个 充 
分 条 件 ， 而 不 是 必要 条 件 。 映 射 四 可 能 不 是 一 种 压缩 ， 但 选 
代 过 程 仍然 是 收敛 的 。 例如， 由 函数 g(x, gy)= 1 十 2y， 


Wz，9) 一 3 十 各 定义 的 映射 @ 不 是 一 种 压缩 ， 因 为 如果 取 


4(8，0) 和 B(8，4) 两 点 ， 我 们 得 到 

7(A, B)=4, ®(A)=(1, 4), ®(B)=(9, 4) 
和 

7( (4)，@(B) ) = 8>1(4, B) 
然而 ， 无 论 我 们 选取 怎样 的 点 M。， 得 出 的 点 列 M。， 如 op 


i … 将 收 1, 4l 
，RMa，… 将 收 伊 于 点 N( 9 可 ，4 襄 )。 


有 时 候 ， 如 果 修 正平 面 上 两 点 间距 离 的 定义 ， 就 能 确定 
选 代 过 程 的 收敛 性 。 实 际 上 ， 两 点 间 的 距离 可 以 有 不 同 的 定 
义 ， 对 于 旅行 者 来 说 ，A、B 两 点 之 间 的 距离 通常 是 计算 他 
从 4 到达 B 所 花 的 时 间 。 图 224 描 述 的 情况 说 明 ，A4、B 两 点 
阅 的 距离 等 于 线段 44、CD 和 DB 长 度 之 和 ( 从 4 点 到 B 点 ， 
必须 先 到 达 桥 CD， 过 桥 后 再 从 D 到 B ) 。 在 一 个 平面 上 ， 如 
果 只 能 沿 着 相互 垂直 的 方向 移动 ( 图 225 ) ， 则 4、B 两 点 间 
的 距离 应 该 定义 为 线段 4C 和 CB 的 长 度 之 和 。 在 另 一 些 情况 
下 ，A、B 两 点 间 的 “距离 ”定义 为 线段 4C 和 CB 中 的 长 者 
是 更 方便 的 。 平 面 上 两 点 间 的 “距离 ”还 可 建立 另外 的 定义 
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( 关于 距离 的 不 同 定义 从 Zu.4. 斯 里 德 (Yu.A.Sckreider) . 
的 《 什么 是 距离 》 一 书 中 可 得 到 详细 的 说 明 ) 。 

两 点 间 的 距离 (4，B) 通 常 具有 下 列 性 质 : 

@@ 任 意 两 点 4、B 间 的 距离 (4，3) 之 0 当 且 仅 当 ;4、 渴 
两 点 重合 时 ，，(4，F)= 0。 

@ 任 意 两 点 4、.5 间 的 距离 满足 对 称 性 : 了 :和 

| 1(A, 5)=7(B, A) a 

@ 对 任意 4、B、C 三 点 ， 满 足 不 等 式 

7(4， Cs<?(4， 3) 十 (CC， 万 ) 

.如果 时 些 对 象 的 集 ， 其 距 高 的 定义 具有 上 述 性 质 ， 出 称 
此 集 为 一 个 度量 空间 ， 它 的 元 案 称 为 这 个 空间 中 的 点 ,度量 
空间 中 的 点 蕉 至 可 能 是 函数 。 对 于 区 间 [ a，& ] 上 的 连 映 站 
数 p(x) 和 和 W(x) 它们 之 苛 的 距离 可 定义 为 |9(*) 一 1 在 
该 区 间 上 的 最 大 值 : 

(9, Y= ,| (3) 一 wz) 

我 们 已 经 看 到 ， 前 面 所 有 定义 的 距离 都 是 满足 条 件 @~ 
@ 的 ， 即 把 平面 转化 成 度量 空间 有 许多 不 同 的 方法 。 
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这 里 有 一 个 非常 有 趣 的 例子 ， 它 定义 的 距离 满足 条 件 @ 
和 @， 但 不 满足 条 件 @@。 假 设 在 山坡 地 带 ，A、B 两 点 间 的 
滤 离 是 由 从 4 走 到 B 的 时 间 来 度量 ， 因 为 上 山 和 下 山 所 花 的 
时 间 不 一 样 ， 所 以 ，71( 4，B) 志 1(B，A)。 
对 于 逐次 逼近 法 解 方程 组 
人 2) 
y=W(x, 5) 
满足 下 面 的 条 件 就 收敛 ， 
上 映射 中 是 由 函数 9(*，2) 和 内 zx，2) 定 义 ， 它 把 区 域 D 映 
射 到 它 自己 的 子 区 域 ， 并 且 至 少 对 某 一 种 “距离 1(A，B) 
是 一 种 压缩 ， 也 就 是 说 ， 存 在 数 4( 0 4 过 1 )， 使 得 对 于 了 
上 任意 两 点 M ,和 M, 满 足 不 等 式 
7( (MI), BM,)) <ar(M1, M,) 
例如 ， 取 函 效 p(x，3)= 1 十 24， 怀 x#，4) 二 3 十 各 ,如 
果 两 点 4(x1，; 1) 布 (x 32) 河 的 距离 由 下 式 定义 : 
C4) (se)+ 2 (ys—y1) 


(99) 


二 当 (x2 一 *1) 一 2 (y,—y 中 


则 由 q(x，2) 和 #4，5) 定 义 的 旬 答 好 是 具有 去 因子 的 压缩 。 
比如 ， 对 于 4(8，0 ) 和 BC 8，4)， 言 1( 有 4， B)=16, 
而 对 于 它们 的 饼 !(4) 和 和 ?(B) 有 
?( Pi), ®(8)) 一 8 


人 么 用 迭代 法 解 方 泣 组 
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3X 一 工 十 23 
3 
y=8+8 


时 ， 虽 然 对 与 之 相关 的 通常 距离 映射 不 是 一 种 压缩 ， 但 迭 
代 过 程 仍然 是 收敛 的 道理 。 


26 .逐次 有 逼近 法 解 线性 方程 
组 的 收敛 检验 


现在 我 们 把 前 面 论述 过 的 收敛 检验 应 用 于 线性 方程 组 。 
通过 选择 不 同形 式 的 “距离 > ， 我 们 将 得 到 各 种 收 伍 检 
验 ， 并 把 这 些 检 验 用 方程 组 的 系数 来 表示 。 
首先 ， 考 虑 两 个 未 知 量 的 线性 方程 组 
{aiix+ais y=b1 2 
(100).. 
人 az :ixX 十 aas3 一 ba: 
设 olis 0，ax: 失 0， 从 第 一 个 方程 中 解 出 x， 第 二 个 
方程 中 解 出 !， 我 们 得 到 一 个 方程 组 


s 
A 


{ QI11 QI11 
y= ba Lay 
G12 A422 
bi 
Qil 


为 简单 起 见 ， 令 


一 Bi， 一 =a1, -6 一 B，， 


om Er 一 Qa， 则 方程 组 为 
人 


人 3 一 az2z 十 B: 
这 样 ， 映 射 中 由 函数 g ( 21 ) 一 xi3 十 B: 和 由 (xs 》 一 Gay 
.101 


(100) 


十 B: 定 义 。 . 
现在 ， 我 们 来 求 系数 a:!、a: 应 满足 什么 条 件 时 ， 映 射 中 
才 是 一 种 压缩 。 
大 家 知道 ,， :4 (zly91 ) 、B(Cxzaygs) 两 点 之 间 的 距离 是 
Y(AB) = (x ox1) :+g —y1) 
且 映 射 中 把 点 4 映射 到 点 41(a19 iB4，az%1 十 B， ) ， 把 点 
了 映射 到 点 B: ( a1y; 十 Bi， Qaxs +B2 )，A1、Bi 两 点 之 间 
的 距离 是 
Y(A1, Bi)= Va ay ) 二 (天 二 x 
一 Wai (ga 一 i) 十 ax 5 一 
用 4 表示 le、ia:| 中 的 最 大 数 ， 即 ， 
4 一 max (ia 中 lazj ) 
则 由 《 101 ) 式 得 不 等 式 
Y(4i 了 ) 和 dz 一 Zi1) 二 0 
=d4Y(4，B) 
明显 地 ， 如 果 <4 1 ， 则 映射 中 是 整个 平面 上 的 一 种 压 
缩 。 倘 若 这 样 ， 正 如 我 们 所 知 逐 次 允 近 法 收敛 。 
于 是 ， 我 们 证 明了 ， 如 果 
max (lx lasl) <1 (102) 
则 逐次 逼近 法 解 方 程 组 ( 100 ) 总 是 收敛 的 。 
现在 ， 我 们 再 把 上 面倒 述 的 收敛 作 一 检验 ， 直 接 用 方程 
组 (100” ) 的 各 项 系数 表示 出 来 。 回 忆 等 式 


a 
Qi 二 一 一 + 
Ql 
a 
aa 一 一 -221 
Aa:2 


将 它们 代入 条 件 ( 102 ) 式 ， 我 们 得 到 下 面 结论 : 
102 


几 逐 次 逼近 法 解 线性 方程 组 
fe 1X 十 al23 一 5 


人 az :xz 十 ozz3 一 Da 
收敛 的 充分 条 件 是 不 等 式 
a :| - 
, 4) <1 


max ( 
满足 . 


这 个 条 件 说 明 ， 方 程 组 各 行 的 系数 中 ， 主 对 角 线 上 的 系 
数 应 该 大 于 非 主 对 角 线 上 的 系数 .例如 解 方程 组 
fx 一 34 一 一 11。 
， 人 6x 十 3 一 10 
了 时， 由 于 这 种 原因 ， 我 们 应 该 从 第 一 个 方程 中 解 出 y， 第 二 
个 方程 中 解 出 zx。 即 


Q12 
a 


有 时 候 ， 六 用 另外 的 与 *、j 成 比例 的 

未 知 量 代 换 方 程 中 的 未 知 量 x 和 y 是 很 有 用 的 。 例 如 考虑 方程 
六 . 
{12x+y=14 

{ (103) 

人 3x 一 24 一 一 1 


对 此 方程 组 


| 
la 

max (2 于 
lei 


a) =max (ss 3) 


所 


3 
2 
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当然 不 满足 逐次 逼近 法 收敛 的 充分 条 件 ， 但 如 果 令 x 一 二 zr 
则 我 们 得 到 方程 组 


Co (104) 
z 一 2g 一 一 】 . 
对 此 方程 组 
ax 人 (全 中 的 引 ) =mox (二 二) 


一 工 。 


这 说 明 ， 方 程 组 ( 104 ) 能 够 用 逐次 百 近 法 求解 。 
当然 ， 解 ( 104 ) 那样 简单 的 方程 组 ， 用 不 着 使 用 逐次 : 
通 近 法 。 但 是 对 于 多 个 未 知 量 的 方程 组 ， 这 种 方法 非常 有 
用 。 我 们 解 方程 组 
/Ci11X1 十 Ci2Xa 十 %。" 十 Ginxn 一 D 
| Qs 1X1 二 as2X2s 十 "十 a2nXn 二 Db 
je 
(anizi 十 anaza 十 .十 annxn 一 pn (105) 
时 ， 收 各 的 充分 条 件 几 乎 与 两 个 未 知 量 的 方程 组 一 样 ， 实 际 : 
上 有 下 面 的 论述 : 
如 果 方 程 组 ( 105 ) 的 系数 满足 下 列 条 件 之 一 ， 则 用 逐 
次 逼近 法 解 方程 组 ( 105 ) 收敛 。 


Dmax laud 


fl ed 1 〈106 ) 
| 
A 
@marD el<1 s (107) 


104 


3 


< (108) 


Ek 


Q@max 3 ai 
j= 全 
次 多 近 法 解 方程 组 ( 105 ) 收敛 。 
在 《106 ) 和 (107 )- 式 中 ， 符 号 《12 的 意思 是 应 除去 
-和 式 中 i=7 的 ' 项 ，( 108 ) 式 中 的 符号 “( 有 ”表示 应 除 盐 
:和 式 中 i=7 和 i=k 的 项 。 要 注意 的 是 ， 如 果 


n/ jn t 
2 到 

‘ 1 

s 7 art 


<1 (109) 
( 其 中 “1” 仍 表示 除去 i=] 的 项 ) 则 条 件 @ 一 定 满足 ， 在 
大 多 数 计算 数学 书 中 ， 条 件 @@ 就 是 ( 109 ) 的 形式 。 
与 前 面 论述 的 一 样 ， 如 果 我 们 要 求 映射 


b a a 
i 1 一 > 一 一 一 二 一 2X。 
QI11 Qi Gil 
b a a 
x np, 
G22 A422 G22 
| b a 
bs 一 bn analzi 一 .， nn 
Qnn Qnn a 


相对 于 某 种 距离 来 说 是 一 种 压缩 ， 则 马上 会 得 到 条 件 ( 106) 
一 (108 )。 洋 际 上 ， 条 件 (106 ) 相当 于 4 (#1，x25 es 
Xn ) 和 B (yi1, 32，… Yn ) 两 点 间 的 距离 为 

TY(4、 了 3 ) 王 max (jz 一 9 和 lz 一 gz …， jzn 一 go| ) 
条 件 (107 ) 相当 于 4、B 两 点 的 距离 为 : | 


Y(4,B)=Z [xyd 


条 件 (108 ) 相当 于 4、 两 点 的 距离 为 “2 
E05 


Y(4BE) = /EC ) 


与 两 个 变量 的 情况 一 样 ， 有 时 用 与 其 成 比例 的 新 的 变量 
代 换 原 方程 中 的 未 知 量 z4，*s，…，xs 是 很 有 用 的 , 即 令 
i=PiXs> y2—=Pp2X2, Yn= pnyny i 
‘pn>0。 
在 这 种 代 换 下 ， 条 件 (106 ) ，( 107 ) ，( 108 ) 分 别 
为 : 


(7 max 并 “laulprc 1 (106’ ) 
i jellatitDhs 本 
(2 ) max 忆 laupsc (107’) 
=atdps 
(3) max eu pe (108 
tl art ps* 


特别 地 ， 当 1 二 la 时 ， 这 些 条 件 成 为 ， 


n/| 
(1) maxD < 1 ( 106'’ ) 
(2”) maxF a - (107’’ ). 
了 {la 
K9l 人 
(38) max 工 2 < 1 (108’’) 
天 ¢9 j= 1077: 


作为 例子 ， 我 们 考虑 方程 组 
: X 一 0.63 一 0.5: 一 一 2.6 
\ 一 0.2x 十 3 一 0.4: 一 3 
\ 一 0.1x 十 0.5g 十 z 一 
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~ 


对 这 个 方程 组 ， 条 件 D、@ 以 及 条 件 (1/)、( 2') 
都 不 满足 。 同 样 地 ，( 109 ) 式 世 不 成 立 ， 因 为 非 主 对 角 线 
元 素 的 平方 和 等 于 1.07， 但 是 


mex 3 kel- 2 6 十 0.5: 二 021 十 0.4 


0.62 十 0.5: 十 0.12 十 0.5:，0.23 十 
0.42 十 0.12 十 0.5: ) 一 0.87<< 上 
所 以 ， 这 个 方程 组 可 以 用 逐次 逼近 法 求解 。 

值得 注意 的 是 上 面 论述 的 所 有 条 件 ， 都 仅仅 是 逐次 逼近 
法 收敛 的 充分 条 件 而 不 是 必要 条 件 。 对 两 点 间 的 焉 离 和 压缩 
条 件 选 择 不 同 的 定义 ， 我 们 就 能 得 出 新 的 收 化 条件， 然而 ， 
在 这 里 我 们 不 继续 往 下 讨论 这 个 问题 了 。 

第 5 节 中 儿 述 过 的 要 点 也 同样 适用 于 线性 方程 组 ， 例 姥 
近似 值 的 最 终结 果 并 不 依赖 于 初始 近似 值 的 选择 ， 所以， 在 
某 次 计算 中 发 生 了 错误 不 会 使 后 面 的 计算 失效 ， 而 仅仅 是 下 
迟 了 向 最 终结 果 逼 近 的 时 间 . 

解 线性 方程 组 ， 逐 次 通过 法 的 形式 是 多 种 多 样 的 . 例如 
在 有 的 方法 里 ， 求 出 近似 值 z2,'"*'? 后 ， 就 以 它 连同 *s'"?， 
Xicn2 m7 一 起 计算 x2 "+1 xs "+19 求 出 以 后 ， 又 以 
x n+ Ka tl9, Am, Xman? 计算 xs +1 如 此 
等 等 。 逐 次 和 逼 近 法 用 于 解 线性 方程 组 的 情况 ， 在 另外 的 书 里 
有 更 详细 的 叙述 。 
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下 A 3 训 人 人 


27 .几何 学 中 的 逐次 通 近 法 


我 们 已 经 令 述 了 逐次 通 近 法 解 方程 和 方程 组 的 问题 。 逐 
次 逼近 法 还 可 用 于 某 些 几何 况 题 ， 例 如 计算 圆周 长 的 问题 w 
正如 大 家 所 知 ， 计 算 圆 的 周 长 实 际 上 是 先 计 算 团 内 接 正 方形 
的 周 长 ， 再 计算 内 接 正 8 边 形 、 正 16 边 形 的 周 长 ， 妇 此 等 
等 。 计 算出 的 这 些 周 长 值 成 一 数列 ， 它 的 极限 值 就 等 于 圆 的 
周 长 。 在 计算 过 程 中 ， 和 次 让 和 全 合作 下 计生 二 的 
结果 。 这 即 是 采用 了 下 面 的 方法 : 、 

用 4 表示 正 2 " 边 形 的 边 ，P 表 示 正 2 " 边 形 的 局 长 
例如 4: 表 示 正 方形 的 边 . 所 以 4: 一 Kv 2 ，P: 王 4RvE 
假设 我 们 已 经 求 出 了 PE， 则 明显 地 有 

Pn 


An= 


在 几何 学 里 已 经 证 明 ,内 接 正 27? 边 形 的 边 a,» 可 以 通过 内 接 
正 n 边 形 的 边 an 和 圆 的 半径 由 下 式 表示 。 
arn=RV 2—V 4 ai/Ri (110) 
因而 ， 内 接 正 2*+! 边 形 的 边 4n+4 可 以 通过 内 接 正 2” 
边 形 的 边 4: 由 公式 表示 ， 即 
Anri=R/ 2—v 4—AiR: 
因为 


A Et 
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出 此 可 得 RR 
i n+1 
Pr+i= 2 Rf 2 /A= 5 

数列 FE:，Ps，…，Fa，… 趋 向 于 圆周 长 度 ， 即 趋向 于 2rR 
的 值 . 所 以 ， 可 以 把 公式 (110) 看 成 是 借助 于 逐次 逼近 法 求 
2xR 值 的 计算 公式 。 用 这 种 方法 可 以 计算 x 的 值 到 任意 位 小 
数 。 

舅 外 还 有 一 种 计算 x 的 近似 值 的 方法 ， 叫 等 周 长 法 。 
这 种 方法 用 一 个 具有 同样 周 长 值 的 正 2**! 边 形 符 换 正 2* 
边 形 ， 假 定 我 们 用 ,表示 正 2" 边 形 的 边 心 距 ， 用 ?7。 表 示 
2" 边 形 外 接 圆 半径 ， 同 样 ， 用 i,+ 1 表示 正 2*+! 边 形 的 边 心 
距 ，7a+1 类 示 正 2 "+! 边 形 的 外 接 贺 半 径 。 

设 4B8 是 内 接 于 半径 为 ?的 圆 的 正 2 * 边 形 的 边 (图 23)。 
将 弧 A 有 8 的 中 点 C 分 别 与 4、B 连 接 ， 并 且 过 三 角形 4CB 的 4C . 
边 和 EC 边 的 中 点 D、E 作 线段 DE。 明 显 地 ， 角 DoH 等 于 能 
Ao5 的 二 六 。 因 此 ，D 上 是 内 接 于 半径 为 oD 的 贺 的 正 2"! 边 


形 的 一 条 边 。 因 为 DE= 寺 4B， 所 以 正 2 "+: 边 形 的 周 长 等 


Pi / 
各 (110’) 


这 个 公式 很 容易 用 三 角 学 导 出 。 明 显 
地 ， 恕 果 a4 是 内 接 正 a 边 形 的 边 ，asn 是 
内 接 正 2 n 边 形 的 边 ( 如 下 图 )， 则 


Pe 
an=2 RsinT,azn=2 Rsin®. 
n 2 


4 ein OQ 1-coso A 
因为 sin 人 -/ 3 ， 内 此 得 由 


一 一 一 二 
aan -2RsinE =2RA/ Leos 了 
2 


从 


vere ey 


图 23 


于 正 2 , 边 形 的 周 长 。 这 即 是 说 ，rny: 一 OD,ls+i 一 DR 
可 以 很 容易 地 计算 出 


1 一 OK 一 二 (111) 
然后 ， 由 直角 三 角形 ODC 又 可 求 出 ，: 
7 一 也 Za (112).. 


公式 CC 111 ) 和 ( 112 ) 就 是 以 fw、 .来 表示 nu 和 ut 的 。 
随 着 z 的 增加 ， 多 边 形 的 周 长 并 不 改变 ， 并 且 数 列 !, 和 
,站内 于 同一 横 限 估 。 _ 这 极限 值 就 等 于 圆 的 半径 ， 该 圆 的 周 


长 就 等 于 多 边 形 的 周 长 。 A 2 
的 多 边 形 ， 则 ?4 和 i 表 向 于 数 二 ， 即 


limy ,一 荆 ， lim¥7 二 一 
了 ~ co 一 > co 人 


例如 ， 如 果 开始 选择 一 个 边 长 为 二 的 正方 形 . 我 们 有 


"= 1 二 . 所 以 ;有 下 面 的 叙述 ， 
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如 果 先 取 ?s 一 2 n=4， 并 由 公式 (111) 和 (112 )》 
计算 715+ :和 1 (nm 一 2，3 ，… ) 。 则 我 们 将 得 出 


lim?， limln_ 
有 一 co n~>o0 x 


这 个 公式 可 以 用 来 求 元 的 近似 值 。 为 计算 出 元 的 值 ， 一 
般 要 渗 续 计算 直到 Y: 和 和 i. 的 值 部 在 要 求 的 精度 之 内 . 此 时 Ys 
或 ,就 是 革 在 精度 要 求 内 的 值 。 


了 1 且 


本 书 使 我 们 了 解 了 逐次 逼近 法 对 制定 计划 、 求 方 根 、 解 
方程 、 计 算 圆周 长 等 各 种 各 样 问题 的 应 用 。 但 所 有 这 些 根本 
不 能 包括 完 逐 次 逼近 法 的 应 用 。 大 量 的 实际 问题 ， 导 出 的 是 
微分 方程 ( 包含 有 未 知 函数 的 微 商 ) 、 积 分 方程 和 各 类 更 复 
杂 的 方程 。 求 解 这 些 方程 的 近似 解 ， 最 有 效 的 方法 之 一 就 是 
应 用 逐次 逼近 法 ， 当 然 ， 就 这 些 应 用 来 说 ， 比 代数 方程 要 复 
条 得 多 。 然 而 ， 可 以 这 样 认 为 ,如 果 不 用 逐次 逼近 法 ,当今 需 
要 着 手 处 理 的 大 量 的 物理 上 和 技术 上 的 问题 ， 没 有 一 个 能 解 
决 。 例 如 ， 计 算 人 造 卫星 运行 轨道 ， 设 计 一 座 原子 能 反应 
堆 ， 研 究 原子 的 结构 等 等 。 不 过 ， 逐 次 逼近 法 对 这 些 问题 的 
应 用 已 超出 了 初等 数学 的 领域 ， 远 远 超过 了 本 书 论述 的 范 
国 。 
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练 习 


下 面 ， 我 们 列 出 一 些 求 方程 近似 解 的 习题 ， 便 于 读者 能 
够 检查 自己 对 书 中 讨论 的 各 种 方法 的 掌握 程度 。 
一 、 用 适 代 法 解 下 列 方 程 *。 


1 42 7 


OO oA Ww nb 
心 
[ 
六 
| 
et 
吕 
3 


© 

| 
3 
wm 
* 


11. x*= 


* 有些 习题 ， 必 须 先 把 方程 改写 为 X=qp《 X ) 的 形式 。 
3 


ul 


14 


13. xz 一 土 V1og(〔(2 十 2 ) 

14. 2%2 一 In (2 十 1 ) 

15. lnx=4—x? 

16. lnx=2—x 

17. x:=ez+2 a 
18. lgx=0.1x 

19, tgx=1gx 


= 1,-z 
20. xz 一 10 


用 牛顿 法 解 下 列 方程 。 
21. xs? 一 5 十 1 一 0 

22. xs 一 9x2 十 20x 一 11= 王 0 

23. 533 一 3x2 一 3x 十 11 一 0 

24. 2%5 十 5x 十 1= 0 

25. sinx 十 x 二 1 

26. x:—1l0lgx—3 二 0 

用 逐次 逼近 法 解 下 列 方程 组 。 要 求 精度 为 0.001. 


sn 0.24—0.1z+0.898 
3 一 0.3x 十 0.15z 十 1.383 
le 一 0.25x 一 0。43 十 3.677 


ww |* A 336 


1 一 一 sin (x—y) 十 0.362 


(ce) 人 0.710 
3 一 1 一 Z 十 1 


练习 答案 


，. 解 : 设 p(x) 一 Er EL (x)= ~ 


2 
《3 和 Gt) 
我 们 有 9 (0) =1>0， 9 (1)= 了 <1， 故 在 区 间 (0， 


1 里 包含 有 方程 的 一 个 根 。 但 是 在 此 区 间 上 不 能 应 用 逐次 遥 
近 法 ， 这 是 因为 | p“(0) | =2>1。 为 此 缩小 区 间 范 围 。 注 


意 到 9 ( 0. 4)=7> 0.4， 所 以 ,方程 的 根 位 于 区 间 (0.4， 
1) 上。 
如 果 0.4 志 x 二 1， 则 |p (*) | si 2<1。 在 此 区 间 


上 可 以 应 用 逐次 允 近 法 ， 取 x 二 0.4， 经 过 11 次 逼近 计 算 后 
得 到 
Xi (xi ) TO0.4655 

所 以 ， 方 程 精度 为 0.001 的 根 是 xz 一 0.4655. 

2. 解 : 设 p (xz) = 二 (x 十 1)5， 风 g(x) 二 3 (x 十 1)*， 
且 9 (一 2 ) = 一 1> 一 2, 9 (一 3 ) = 一 8 之 一 3。 所 以 在 区 
间 ( 一 83， 一 2)] 上 有 方程 的 根 ， 然 而 在 此 区 间 上 ,| (*) | 
> 1， 不 能 应 用 逐次 百 近 法 。 

把 方程 改写 为 

x 一 以 x 一 1 


卫 15 


刻 有 友 〈(x ) 一 吕 2 一 1， 由 (2 ) 一 3 在 区 间 ( 一 3 


—2) 上 , |Y < 了 可 以 用 逐次 逼近 法 。 ; 


. 取 x, 二 一 2， 得 到 ze 一 和 (xze) ~ 一 2.325。 
所 以 在 精度 0.001 内 ， 方 程 的 根 为 x 二 一 2.325。 


六 


3. 解 : 设 p (xz ) 一 4 十 人 ， 则 


q(x*)= 


i +1)’ 
图 24 表 明 ， 直 线 y 一 x 与 曲线 9=4 十 * /< 二 有 两 个 实 


点 ， 分 别 在 区 间 ( 一 1，0) 和 (4，5) 上 。 下 面 逐一 进行 如 
- 解 。 


在 区 间 (4，5] 上 ， 有 19“ (2 ) | 56< 1 取 x， 


二 4， 计 算得 *s 二 p (xs ) 久 4.870。 所 以 在 精度 0.001 内 有 


| 


-十 
1 
! 
1 
1 


16 


滩 . 


$b 


MR A DE A [= 


X=4.870。 
在 区 间 [ 一 1，0)] 上 ， 不 能 直接 用 逐次 逼近 法 ， 把 方程 
正 写 为 - 


ys 2 一 上 -: 
| (x 4) zi 
由 此 得 
UR 
二 [= 二 x 十 1 
1 
即 “TT 1 


这 里 ， WT 1 (4)= 1 性 去 . 明 地 ， 


对 一 1 Kx<0 有 |W(*) 1 < 直 <1 可 以 用 梁 次 


逼近 法 求解 。 取 x, = 二 0， 得 *, 二 (x ) 久 一 0.9840， 即 在 糖 
度 0.0001 内 ， 有 2Y 一 一 0.9840。 

所 以 ， 我们 求 得 方程 的 两 个 根 是 += 一 0. 9840 和 x 二 
4.870。 LA 

4- 解 设 p1 (x*) 一 2 十 x ,qi 二 2 一 xw。 由 图 
25 可 以 署 出 ; r= 2 + 在 区 间 (3，4) 上 有 一 根 。 
在 此 区 间 上 : 


lewW | -<! a 


可 以 用 逐次 逼近 法 。 取 zi 一 4， 
得 x+=qi (2 ) 尽 3.353， 因 
此 ， 在 精度 0.001 内 x 二 3.353。 

现在 解 方程 * 一 2 一 x 。 
它 的 根 显 然 是 xz 一 1 。 图 25 
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所 以 方程 的 根 是 zx 一 3.353 和 xz 一 
5. 解 ; 这 里 ，9(z) 一 以 5 二 29/(z) 一 


且 有 9(1)=&4<2，9(2)=&3>1。 
所 以 ,方程 在 区 间 (1，2) 上 有 一 根 。 在 此 区 间 上 


19’ (x*) | < <ava<l. 取 xzi: 一 1， 计算 得 zs 一 9(ze》 
X21.516。 

因此 ， 在 精度 0.001 内 ， 方 程 的 根 是 xz 一 1.516。 

6 ， 解 ， 把 方程 改写 为 

X=arctg ( 4 一 zx ) ， 5 

在 这 里 , p (*) 一 arctg ( 4 一 *), 9(1)=arctg3~1.25r 
9 (2)=artctg2~1.10。 可 见 ， 方程 在 区 间 (1， 2) 上 有 一 
根 。 


在 此 区 间 上 ,| 9’ (x) | 一 


rs 
3 8/(5—x):, 


次 逼近 法 。 取 *, 二 1 ， 得 *, 二 9 (*。 ) 沪 1.225。- 

所 以 ， 在 精度 0.001 内 ，x 一 1.225。 去 

了 7. 解 ， 方程 有 一 根 为 x= 0 。 

从 图 26 可 以 看 出 ， 方 程 的 另 一 根 是 正 数 ， 所 以 应 满足 方 
程 交 
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2X 一 wsinx 


一 vsinz ‘ 一 CO 因 闷 
这 里 ,9(2) 二 wsinz， q (0) =n 因为 


9( 证 ) 一 Vsia 二 <Y5.5764> 评 
9 (1)=VsinlYV0.8414< 1 . 
所 以 ， 方 程 在 区 间 [ 去 ,1 ) 上 有 一 根 . 在 此 区 间 上 


1 
Cos 
19'() |< 2 ~0. .8703 一 1 


2 fanl 1.3846 
Sn2 


逐次 逼近 法 收敛 。 取 xi 二 1， 得 到 x*1 = 二 p(x1) 二 0.8768。 所 
以 ,在 精度 0.0001 内 ,方程 的 第 二 个 极 是 zx 一 0.8768。 

8 和. 解 ， 这 个 方程 的 解法 与 上 题 一 样 。 先 把 方程 写 为 

- X 一 只 sinxX 

再 取 x, = 1 ， 计 算得 到 ze=9 ( x*。) 六 0.9286。 房 疝 方程 
的 一 个 根 等 于 0.9286， 精 度 为 0.001。 

因为 方程 两 边 都 是 奇 函数 ， te 
一 0.9286。 

方程 的 第 三 个 根 是 xz= 0 。 

9. 解 ， 原 方程 可 以 写 为 


-一 sinx (-<*<#) 
从 图 27 可 以 看 出 ， 方 程 在 0 入 之 间 有 一 单 根 。- 


1 1 


这 里 ，qg(x) 二 arcsin” GT 


2 9 (x*) 一 
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在 区 间 (0， 已 ] 上 有 
1 
A 
lo’) 1 < a<1 


取 x, 二 0， 得 *, 二 p(x， ) 多 0.3422。 

所 以 ，x 一 0.3422 是 方程 具有 精度 0.0001 的 根 。 

10， 解 ;因为 cos0 = 二 1，cosi 之 0 .所 以 ， 方向 00 时 
在 区 间 (0，1) 上 有 一 根 。 

在 此 区 间 上 ， | 9g'(x) | 和 sin1< 1 ， 可 以 用 逐次 逼近 
法 。 取 x, = 1 。 计 算得 到 x 二 0.7391， 具 有 精度 0.0001。 

11. 解 ， 从 图 28 可 以 看 出 ,方程 的 正 根 位 于 函数 3 一 cosx 


的 图 像 与 4 轴 的 交点 近 旁 ， 且 离 它 左 边 的 交点 是 瑟 十 24m， 
右边 的 交点 是 全 十 ( 2k 十 1)x。 为 了 求 出 方程 在 + 一 后 十 xz 
邻近 的 根 ， 令 x 一 ax 一斑 一 纺 则 方程 为 


4 一 一 一 一 一 一 (人 一 1 二 


cos(y 十 nx 十 分 ) ai 
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图 28 


Jy， 所 以 ， 方程 又 可 写 为 


==(— 多 机 了 1 Cs 
y=( 1) yi arcsin GFTnxt 到 
这 样 . 
gg(y)=(~—1 ) +iaresjin 1 = 
.十 NX 十 pm 


Ch) ". 


9' (4 ) = 一 一 一 和 Re 
Cot 玫 tm 一 1(gtnxt 茹 ) 
显然 ， 在 点 I= 0 附近 有 a 
lq9’(9) |<a<1: 
可 以 用 逐次 逼近 法 。 为 了 求 出 n== 工时， 精度 为 0.001 的 根 ， 
取 y。 二 0， 计 算得 y, sg (ys ) 一 0.204 
所 以 ，y~0.204。 由 此 得 x 二 二 44.917 


为 子 求 出 第 一 个 负 根 ， 取 n= 一 1 ,- 我 们 得 到 方程 
121 


1 


和 


ge 一 9 (ye ) 祥 一 0.503. 

所 以 ，y 污 一 0.503， 由 此 得 * 污 一 2.074 

对 | | 的 较 大 值 ， 逐 次 百 近 法 给 出 了 一 个 计算 y 值 的 近 
似 公式 ， 


y=arcsin 


2 取 y。 二 0， 计算 得 


1 
yo(y0)=(—1)"tarcsin 


x 
nx 二 二 
E34 Pr n+l 疙 i 
(一 1) (2n+1)7x 
，X 包工 。 ee ER ” 
xs (2r+1) 二 (一 1) CT 元- 


12， 解 : 取 x: 一 0， 计算 得 到 xs 一 中 (xs ) 1.088c 
所 以 ， 方程 具有 精度 0.001 的 根 是 x 二 1.088。 
13， 解 ， 先 解 第 一 个 方程 


z=,/1g (x+2) 
lge 

2 (x+2(VIg (x+2) 
因为 9 (0) = 二 Vig2>>0, 9 (1) = 二 Vlig3<<1。 故 方程 在 区 
间 (0y 世上 有 一 根 。 在 此 区 间 上 ，| 9 (zx) | <&<1, 因 
此 能 用 逐次 逼近 法 。 取 x, = 1 ， 得 xs ~ 《xs ) 0.6507。 
所 以 , 方程 z=Wig tx 二 2) 具 有 精度 0.0001 的 根 等 于 
0.6507。 

再 解 第 二 个 方程 
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这 ,CD IEG 本 ,01 (x)= 


3X 一 一 lg(x 十 2 ) 

在 这 里 ，9(2) 一 tv lg(z+2), «iB 55， 
9( - 划 = -vi 谭 5= 一 0.42" 芍 劣 程 在 区 间 ( 去 ，0) 上 
有 二 根 。 取 x 二 0， 计算 得 到 zs 二 9 (x4 ) 心 20.4397。 所 
以 ， 在 精度 0.001 内 有 x 一 一 0.4397。 1 

14. 解 :方程 有 一 个 根 是 *= 0 。 为 求 其 余 的 根 ， 把 方 
程 写 为 

X 一 土 MIln(x 十 1 ) 


对 于 方程 <=wVin (x 二 TI )， 有 
$$)= /In (这 ) > 二 (1)=vIE5<1 
因此 ， 方 程 在 区 间 ( 寺 ， 1) 上 有 一 根 。 凤 为 
/ 人 Se 
9 《455TTJVETST 
所 以 在 区 间 ( 却 ，1) 上 ，19'(z) 1 <a< 1。 取 zi=1， 
得 加 =9 (xx ) 0.7469, 所 以 ， 3 一 0， 7469 具 有 精度 0. 0001, 


对 于 方程 zx 一 一 win (x 十 1 )， 除 z= .0 外 没有 另 外 的 
根 。 

所 以 ， 方程 的 要 是 x 二 0 或 + 二 0.7469。 

15. 解 ， 把 方程 改写 为 

2X 一 4 一 lnx 
这 里 9 (4) = 一 ia 一 inaz， 9 9 (1)=2, (2) 一 4 二 1n2， 
一 1 
VDT 本 全 
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图 29 图 30 


因为 p (1 ) =2，9 (2 ) =vV4 一 In2 之 2; 故 方程 在 区 间 (1， 
2 上 有 一 根 。 从 图 29 中 明显 知道 ， 方程 没有 另外 的 根 。 取 zf 
二 2， 我 们 得 到 x 入 9 (2x4 ) N1.841。- 

所 以 ， 具 有 0.001 的 根 是 2 一 1.841。 

16. 解 ， 把 方程 写 为 


X 一 2 一 Inx 


的 形式 。 这 里 9 (4) 一 2 一 Inzr 9/(z ) 二 一 二 ;1 从 图 30 可 


bd 
以 看 出 坊 程 的 根 在 区 间 (1;2) 上 ,在 此 区 间 . 上 ， | 9 (3z)| 委 : 
1。 取 x1 二 1. 5， 我 们 得 到 xssg (xis ) 久 1.557。 
所 以 ，x==1.557 具 有 精度 0 .001。 
17， 解 ， 从 图 31 中 可 以 看 出 ， 方 程 仅 有 一 负 根 。 把 方 禄 
写 为 : 


. 
* 生 一 Ver 十 2 


则 -9(>) 一 一 ez 十 3，9/ 6 Ey 六 ， 
9(—1)=—vVe-!+2~—1.54, 9(—2)=—VE7T 
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46 - 
所 以 方程 的 根 位 于 区 间 ( 一 2， 
一 上。 到 x, 一 一 1, 计算 得 
到 Xs (x, ) OT —1.492。 

所 以 ， 在 精度 0.001 内 有 
x 一 一 1.492。 

18. 解 . 明显 地 ， 方 程 有 
一 根 是 xz=10. 为 求 方程 的 第 二 
个 根 把 方程 写 为 :zx 一 102 1!= 图 31 
这 里 ，9 (x*) =10%"!*, qg‘(x)=0.1Xx10*"*rtn19, 且 
9p(1) =10'"! 泡 1，q (2 ) =10" :<2。 因 此 方程 在 区 间 
《1，2) 上 有 一 根 。 在 此 区 间 上 

Il9’ (x) |<0.1xi0°*1n10~0.37<1 

能 用 逐次 逼近 法 。 取 x，, = 二 .2 ， 得 x*1 守 q(x; ) 1.372。 
所 以 ，x 二 1.372 有 精度 0.001。 i 


19. 解 ， 从 图 32 豆 以 看 出 ,方程 在 每 个 区 间 ( 竺 十 nx。 和 


125 


十 (mn 十 1 ) r)] 4 一 0，1，2 。。。 上 都 有 一 根 。 且 且 这 些 根 都 在 
每 个 区 间 的 右 其 。- 为 了 求 第 一 个 正 根 ， 作 代 执 2 一 号 一 六 “ 
则 方程 为 i 
ct8g 一 所 (六 x 一 g ) Us 
由 此 得 。 ，. 
y=arcctg [as ( 3 一 起 ] 《0<3g<z)》 


这 里 ，9 (3 ) 一 arcctg [1 ( 3r—s ) ] 


5 一 1ge 
(19 = {1pie: (Ba) ] (3r—s) 


在 区 间 (0，x) 上 , 方程 有 一 根 , 且 | 9^( 8) [<1. 应 用 逐次 
逼近 法 ， 取 yi 二 0， 得 gsg (y。 ) 1.059。 
所 以 ， 在 精度 0.001 内 ， 有 y= 二 1.059 必 基 x 二 3,654。 


， 为 了 求 第 二 个 正 根 ， 令 x 一 x 一 y， 则 方程 为 


y=arcctg [as ( Sry )] 


取 y = 0， 计算 得 34sg (ys) 久 0.870。 
所 以 ， 在 精度 0.001 内 ， 有 3 一 0.870， 即 xz 一 6.984。、 
20， 解 ， 从 图 83 可 以 看 出 ,方程 在 0 和 1 之 间 有 一 根 。 


这 里 ， 9 (7) = 9’(*) 一 一 溃 c…, 且 在 区 间 (0，1Y; 
上 ， 不等式 1 9’(*) | 生硬 成立， 完全 可 以 用 逐次 通 近 法 
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求解 。 取 2i 一 0， 有 zs (xs ) 必 0.091.. 
所 以 ， 在 精度 0.0001 内 ，x 一 0.091。 
21.. 解 : 今 1 (x ) = 二 x: 一 5x 十 1， 则 

f(x)=3x—5, 1/1’'(x)=6x 

使 用 牛顿 公式 有 | 


a Bi—5Ba,, 
Bari—Bn— Ei 


函数 在 一 些 点 上 的 值 如 下 表 。 


从 表 中 可 以 看 出 , -方程 :一 5x 十 1 在 区 间 ( 一 83， 一 2)， 
(0，1) 和 (2，3) 上 各 有 一 根 。 | 

先 求 区 间 ( 一 3， 一 2) 上 的 根 , 在 此 区 间 上 ， 因 为 
1'' (x*) 过 0， 所 以 ,我们 取 初 始 值 8。 二 一 3( 因为 1(B。》 
-一 11<<0 ) 。 计 算得 
-3 人 ttl. 
继续 计算 得 B, 8B. 盖 一 2.331。 所 以 ， 在 区 间 ( 一 3， 一 2) 
土 ， 方 程 精度 为 0.001 的 根 是 xz 一 一 2.331。 

其 次 求 区 间 [0, 1) 上 的 根 .在 此 区 间 上 ,1'‘(x) 之 0， 
因此 取 B ,= 0。 计算 得 


1 二 


Bs ~B YO0.202 _ i 
所 以 ， 在 区 间 ( 0， 1) 上 ,方程 的 根 是 x 二 0.202， 具 有 精度 
0.001。 
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最 后 求 区 间 (2， 3] 上 的 根 。 取 B。 一 3， 计 算得 
pr 3 409 
继续 计算 得 
BBs OT2.128 5 
所 以 ， 在 区 间 (2，3] 上 ， 方程 的 根 是 * 二 2.128， 精 度 为 
0.001。 
”于 是 我 们 求 出 了 方程 的 三 个 根 ，xz4 一 一 2.331，xa 一 
0.202，xs 一 2.128。 ES 
再 用 改进 的 弦 截 法 解 这 个 方程 。 在 区 阅 ( 一 8.72)t, 
我 们 求 出 和 
i 


dt "F372 


= 3+11= Si 214 


在 此 区 间 上 ， 因 为 1”” (4) <0, 所 以 表 数 间 线 是 向 下 加 
的 ， 才 由 下 式 计算 qt ee 

ai 一 3 一 (一 3 ?i 3 7 —2.214) ~2.298 
接着 再 计算 人 


二 -2.293 十 2.214 
oa 2.293—/(—2.293) gi) /C2 a1 


久 一 2.331 ， 

这 结果 与 上 面 得 到 的 * 值 一 样 ,也 在 精度 0.001 内 。 
使 用 同样 的 方法 还 可 求 出 方程 在 区 间 ( 0， 1 和 (2 s) 
上 的 根 。 


22。 解 ， 这 里 ， 我 们 有 
1 (x*)=x— 9x:+20x—11 
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1’ (x*)=3x:—18x+20 
1''(*)=6*—18= 0 
” 计算 函数 值 如 下 表 : 


f(x) 


-11 


所 以 ， 方程 的 根 在 区 间 (0，1), (2，3) 和 (5,. 6] 上 ， 
在 区 间 (0 , 1] 上 , 取 B。 一 0 ,计算 得 到 BexBis0。 834， 
在 区 间 (2 ，3) 上 ， 取 B。= 3 ， 计 算得 到 B:wB,2.216; 
在 区 间 (5，6)】 上 ， 取 B= 6 ， 计 算得 到 B. 必 Be 冯 5.249。 
所 以 ， 我 们 求 出 方程 的 根 是 ，x: 一 0.834，x: 一 2.216， 
Zs 三 5.249， 精 度 为 0.001。 
23. 解 。 这 里 
/1(z) ee-3z 30+11 
j (zx) 一 3x: 一 6x 一 3 
11/’’(x)=6x—6=6 (x—1) 


计算 函数 值 如 下 表 : 
x -2 | -1 | 。 | 1! | 2 3 
f(x) | ~s wo|u|e| 1 |2 


所 以 ， 方程 在 区 间 ( 一 2， 一 1 ) 上 有 一 实 根 ,为 求 此 根 ， 我 
们 取 B ,二 一 2 ， 计 算得 、 
B: 祥 8 之 一 1， 847 
所 以 ，x 一 一 1.847， 精 度 为 0.001。 
24。 解 ， 此 题 中 : 2 
. 1 (x)=x:+5x+1 ' 
f(x)=52:+5 


jx ) =20x* 
列 函数 值 表 ， ， 
-1 | 0 1 因此 , 方程 在 区 税 ( 一 区 0) “ 
fz)| -5s | 1 | 7 | 上 有 有 一 根 。 取 B 一 一 1 在 精度 


0.0001 内 | 
B:sB. 必 一 0.1999 
所 以 ，x 一 一 0.1999 满 足 精 度 0.0001。 
25. 解 ， 这 里 有 
/(2z ) 一 sinx 十 zx 一 】 
f/f’ (x ) 一 cosz 十 1 
1’’(*)=—sinx 


f(x ) 的 函数 值 如 下 表 ， 8 
3 0 1 2 因此 ， 方程 在 区 间 (0， 1) 
f(x) | 一 1 | 0.8115 | 1.9093 | 上 有 一 根 。 取 B。 == 0 ,我 们 求 
得 B: 多 ps 六 0.5110。 
所 以 ， 在 精度 0.0001 内 有 x 一 0.5110。 和 
26. 解 ; 在 这 里 


t(x ) 一 22 一 10logx 一 3 


/ 10 
?t“(x) 22 一 ni0 _ 


1 (4#) = 2 ++io 
f (2% ) 的 函数 值 表 为 : 
0.5 1| 2 因此 ， 方 程 的 根 位 
o.26 | -2 | -mms 1) 和 (2， 


bE3 


f(x) 
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ss 


3) 上 。 

在 区 间 (0.5， 1) 上 , 取 B。 = 二 0.5， 得 到 B ,Bs 心 0.535. 
所 以 ， 满 足 精度 0.001 的 根 等 于 0.535。 

在 区 间 (2， 3] 上 ， 取 B ,= 二 3， 得 到 B:sBs 关 2.705。 

所 以 ， 方 程 有 两 个 根 ， ,一 0.535，x :一 2.705。 

27。 解 : (a ) 取 xo 二 0， yo 二 0, 20= 0。 经 过 几 次 
逼近 计算 后 ， 得 到 具有 精度 0.001 的 近似 解 是 


YX 一 1.021 
3 一 2.150 
z=3.072 


(b ) 取 x。== 0，yo 二 0。 计算 几 次 以 后 的 近似 解 ( 精 
度 为 6.001 ) 是 : 
让 
y=0.310 
(c) 取 x*,=0, y=0, 计算 得 
人 


3 一 2.000 
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